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Abstract 

Iron  (Fe)  is  an  essential  micronutrient  for  marine  phytoplankton  often  limiting  phy¬ 
toplankton  growth  due  to  its  low  concentration  in  the  ocean  and  thus  playing  a  role 
in  modulating  the  ocean’s  biological  pump.  In  order  to  understand  controls  on  global 
Fe  distribution,  the  decoupling  between  Fe  and  PO4  and  the  sensitivity  of  surface 
nutrient  concentrations  to  changes  in  aeolian  iron  supply,  I  use  a  hierarchy  of  ocean 
circulation  and  biogeochemistry  models. 

I  formulate  a  mechanistic  model  of  iron  cycling  which  includes  scavenging  onto 
sinking  particles  and  complexation  with  an  organic  ligand.  The  iron  cycle  is  coupled  to 
a  model  of  the  phosphorus  cycle.  The  aeolian  source  of  iron  is  prescribed.  This  system 
is  examined  in  the  context  of  a  highly  idealized  box  model.  With  appropriate  choice 
of  parameter  values,  the  model  can  be  brought  into  consistency  with  the  relatively 
sparse  ocean  observations  of  iron  in  the  oceans. 

I  implement  this  biogeochemical  scheme  in  a  coarse  resolution  ocean  general  circu¬ 
lation  model,  guided  by  the  box  model  sensitivity  studies.  This  model  is  also  able  to 
reproduce  the  broad  regional  patterns  of  iron  and  phosphorus.  In  particular,  the  high 
macro-nutrient  concentrations  of  the  Southern  Oceans  result  from  iron  limitation  in 
the  model. 

I  define  a  tracer,  Fe*  that  quantifies  the  degree  to  which  a  water  mass  is  iron 
limited.  Surface  waters  in  high  nutrient,  low  chlorophyll  regions  have  negative  Fe* 
values,  indicating  Fe  limitation,  because  aeolian  surface  dust  flux  is  not  sufficient  to 
compensate  for  the  lack  of  iron  in  upwelled  waters. 

The  oceanic  residence  time  of  Fe  is  ~285  years  in  the  model,  confirming  that 
transport  plays  an  important  role  in  controlling  deep  water  [Fej-].  Globally,  upwelling 
accounts  for  40%  of  ’new’  iron  reaching  the  euphotic  zone. 

Due  to  the  potential  ability  of  iron  to  change  the  efficiency  of  the  carbon  pump 


3 


in  the  remote  Southern  Ocean,  I  study  Southern  Ocean  surface  P04  response  to 
increased  aeolian  dust  flux.  My  box  model  results  suggest  that  a  global  ten  fold 
increase  in  dust  flux  can  support  a  P04  drawdown  of  ~0.25/xM,  while  the  GCM 
results  suggest  a  P04  drawdown  of  0.5  /M. 

Thesis  Supervisor:  John  C.  Marshall 
Title:  Professor 

Thesis  Supervisor:  Edward  A.  Boyle 
Title:  Professor 
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Chapter  1 


Introduction 


Iron  (Fe)  is  known  to  be  an  essential  micronutrient  for  marine  phytoplankton  capa¬ 
ble  of  limiting  phytoplankton  growth  due  to  its  low  concentration  (Coale  et  al,  1996; 
Martin  et  al.,  1994;  Price  et  al,  1994).  Ocean  regions  with  such  low  Fe  concentrations 
are  typically  characterized  by  high  surface  concentrations  of  macronutrients,  but  low 
chlorophyll  (referred  to  as  HNLC  regions).  The  Southern  Ocean,  northern  North  Pa¬ 
cific  and  equatorial  Pacific  are  classified  as  HNLC  regions  (Figure  1-1).  Fertilization 
experiments  have  shown  primary  productivity  in  these  regions  responds  to  Fe  addi¬ 
tions  (Martin  et  al., 1994;  Coale  et  al.,  1996;  Boyd  et  al.,  2000).  At  their  outset  in  the 
North  Atlantic,  the  deep  waters  are  replete  in  iron  relative  to  PO4.  However,  due  to 
the  different  biogeochemistry  and  decoupling  of  P04  and  Fe  in  the  deep  ocean,  when 
the  waters  upwell  again  (e.g.  in  the  Southern  Ocean  and  the  northern  North  Pacific), 
they  are  depleted  in  Fe  relative  to  phosphorus.  Unless  the  external  aeolian  supply  is 
sufficient  to  offset  the  deficit,  iron  is  limiting. 

Evidence  from  ice  cores  (Petit  et  al.,  1999)  and  sediments  (Rea,  1994),  in  addition 
to  suggestions  from  numerical  models  (Mahowald  et  al.,  1999),  indicate  an  increased 
global  aeolian  dust  supply  of  2-5  times  during  the  Last  Glacial  Maximum  (LGM) 
and  an  increase  as  high  as  20  times  at  high  latitudes.  An  increased  dust  flux  might 
have  increased  nutrient  utilization  (Francois  et  al.,  1997;  Kumar  et  al.,  1995)  and 
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Chlorophyll  a  Concentration 
mg/m3 


Figure  1-1:  Top  panel  displays  the  surface  P04  from  Conkright  et  al.  (1994)  in  pM. 
The  bottom  panel  displays  the  yearly  averaged  chlorophyll  a  concentration  (mg/m3) 
as  sensed  by  SeaWiFs  (Figure  courtesy  of  NASA). 
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subsequent  drawdown  of  pC02  in  the  Southern  Ocean  during  the  LGM.  However, 
this  hypothesis  is  yet  to  be  clearly  supported.  While  the  field  experiments  have 
clearly  shown  an  increase  in  primary  production  following  iron  fertilization  in  the 
HNLC  regions  (Martin  et  al., 1994;  Coale  et  al.,  1996;  Boyd  et  al.,  2000)  there  is,  as 
yet,  no  clear  evidence  of  a  change  in  export  production  which  would  be  necessary  for 
modulation  of  the  biological  pump  of  carbon. 

While  iron  may  play  a  significant  role  in  regulating  primary  productivity  and 
biological  drawdown  of  C02,  the  processes  controlling  its  distribution  in  the  global 
ocean  are  poorly  understood.  Yet,  companies  have  formed  (Lam  and  Chisholm,  2002) 
with  the  aim  of  fertilizing  the  ocean  with  iron  to  mitigate  rising  C02  levels  in  the 
atmosphere  due  to  fossil  fuel  burning. 

Much  of  the  emphasis  of  iron  modeling  studies  have  focused  on  the  cycling  of  Fe 
within  the  ecosystem  in  the  upper  ocean  (Armstrong,  1999;  Christian  et  al.,  2002 ; 
Leonard  et  al.,  1999;  Moore  et  al.,  2002).  These  models  must  impose  the  upwelled 
iron  reaching  the  euphotic  zone.  Only  two  studies  have  focused  on  modeling  the  large 
scale  distribution  and  maintenance  of  Fe  in  the  deep  waters  (Archer  and  Johnson, 
2000;  Lefevre  and  Watson,  1999),  the  source  for  upwelled  iron.  By  building  on  these 
previous  deep  water  iron  studies,  incorporating  new  understanding  and  field  data,  this 
thesis  focuses  on  understanding  what  controls  the  oceanic  distribution  and  supply  of 
iron  to  the  euphotic  zone  and  the  effect  increased  dust  flux  has  on  nutrient  drawdown 
in  the  high  latitudes.  I  test  various  parameterizations  of  iron  cycling  based  on  its 
biogeochemical  properties  within  the  context  of  a  simple  six-box  model  and  a  physi¬ 
cally  more  sophisticated  three-dimensional  coarse  resolution  ocean  general  circulation 
model.  I  find  that  the  binding  of  Fe  to  an  organic  ligand  is  able  to  counteract  the 
loss  of  Fe  due  to  scavenging.  The  dominance  of  scavenging  over  transport  at  detph 
leads  to  the  decoupling  between  Fe  and  PO4  in  the  deep  ocean. 
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Modern  Aeolian  Flux  (mg  Fe/m2/yr) 


Figure  1-2:  Average  yearly  Fe  flux  extrapolated  from  in  situ  marine  boundary  layer 
measurements  by  Gao  et  al. (2001). 

1.1  Biogeochemistry  of  Iron  in  the  Oceans 

Like  other  metals,  such  as  lead  and  aluminum,  iron  has  an  episodic  aeolian  source 
to  the  surface  ocean.  Global  estimates  of  annual  atmospheric  iron  deposition  range 
between  14-32  x  1012  g  Fe  (Ducc  and  Tindale,  1991;  Tegen  and  Fung,  1995;  Ma- 
howald  et  al.,  1999;  Gao  et  al.,  2001)  with  deposition  strongest  in  the  North  Atlantic 
and  Indian  Oceans  (Figure  1-2).  Dust  deposition  is  lowest  in  the  South  Pacific  and 
Southern  Ocean. 

Assessing  the  bioavailability  and  solubility  of  aeolian-derived  Fe  is  a  major  research 
focus.  The  extent  of  dissolution  in  seawater  appears  to  depend  on  the  aerosol  source, 
degree  of  atmospheric  processing,  residence  time  of  aerosol  Fe  in  the  euphotic  zone, 
particulate  load  present  in  surface  waters,  and  the  chemical  reactions  aeolian-derived 
Fe  is  subjected  to  in  surface  waters  (Jickells  and  Spokes,  2001).  Recent  studies  suggest 
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the  solubility  of  Fe  in  dust  may  be  below  5%  (Jickells  and  Spokes,  2001;  Spokes  and 
Jickells,  1996). 

Other  recent  evidence  indicates  that  Fe  is  found  in  both  the  soluble  (<  0.02  pm) 
and  colloidal  size  classes  (0.02-0.4  pm)  (Wu  et  al.,  2001).  Previously,  due  to  analytical 
difficulties,  investigators  were  unable  to  distinguish  between  these  two  forms  of  Fe. 
Colloidal  iron  particles  may  be  less  bioavailable  than  soluble  iron  (Hudson  and  Morel, 
1990;  Wells  et  al. ,  1983)  and  could  aggregate  into  larger  particles  that  sink  out  of 
the  water  column  (Honeyman  and  Santschi,  1989).  Since  most  of  the  available  Fe 
measurements  were  not  able  to  distinguish  between  soluble  and  colloidal  Fe,  I  do  not 
consider  the  role  of  colloids.  Rather,  in  this  thesis,  ‘dissolved’  Fe  refers  to  Fe  that  has 
passed  through  a  0.4  pm  filter.  As  we  learn  more  about  the  nature  of  Fe-colloids  in 
the  ocean,  I  hope  to  incorporate  this  knowledge  in  the  future. 

Iron  is  removed  from  the  water  column  by  scavenging  onto  sinking  particles.  Direct 
quantitative  estimates  of  scavenging  rates  of  Fe  have  not  yet  been  made,  though 
Bruland  et  al.  (1994)  indirectly  estimate  a  residence  time  for  Fe  with  respect  to 
scavenging  between  70-140  years  in  the  water  column.  Thorium  (Th)  is  a  metal 
that  has  similar  abiological  properties  to  Fe.  Bacon  and  Anderson  (1982)  calculate 
an  oceanic  scavenging  rate  for  Th  and  also  suggest  that  scavenged  Th  is  released 
back  to  the  water  column.  They  describe  the  latter  process  as  a  first  order  reaction 
proportional  to  the  particulate  Th  concentration,  estimating  redissolution  rates  of 
1.33-6.30  yr_1.  Since  Fe  and  Th  have  similar  metallic  properties,  it  seems  reasonable 
to  speculate  that  scavenged  Fe  on  particles  may  also  be  subject  to  redissolution. 

The  vertical  profile  of  iron  reflects  its  role  in  the  biological  cycle  and  atmospheric 
source.  Its  vertical  profile  falls  into  two  categories  based  on  the  depth  of  the  mixed 
layer,  rate  of  biological  productivity  and  dust  flux.  Regions  with  deep  mixed  layers, 
high  productivity  and/or  weak  dust  flux  are  characterized  by  depleted  [Fer]  at  the 
surface  that  increases  with  depth  (Figure  1-3 A).  In  regions  with  stratified  mixed 
layers,  low  biological  productivity,  and  strong  dust  flux,  iron  builds  up  near  the 
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Figure  1-3:  Vertical  profile  in  the  A)  Equatorial  Pacific  at  3°S,  140°W  (Johnson  et 
al.,  1997)  and  B)  western  North  Atlantic  at  34.8°N,  57.8°W  (Wu  et  al.,  2001). 

surface,  has  a  minimum  near  the  base  of  the  mixed  layer  and  then  increases  with  depth 
(Figure  1-3B).  In  Figure  1-4A,  I  present  a  compilation  of  surface  [Fer]  from  published 
and  unpublished  measurements  courtesy  of  E.  Boyle.  In  high  dust  flux  regions,  such 
as  the  North  Atlantic,  surface  [Fer]  can  be  elevated.  Due  to  the  analytical  difficulty 
of  measuring  iron,  the  deep  water  iron  distribution  is  currently  poorly  resolved,  but 
it  is  clear  that  large  scale,  deep  water  Fe  gradients  do  not  mirror  those  of  nitrate  and 
phosphate.  Rather  concentrations,  at  ~1000m,  are  highest  in  the  Atlantic  (  0.6-0.8 
nM),  intermediate  in  the  Indo-Pacific  basin  (0.4-0. 7  nM),  and  lowest  in  the  Southern 
Ocean  (0.2-0.3  nM)  (Figure  1-4B).  For  depths  greater  than  2500m,  few  measurements 
exist  (Figure  1-4C),  but  suggest  highest  [Fer]  in  the  Atlantic  basin  and  lower  [Fer] 
in  the  Southern  Ocean  and  Pacific  basin. 

Field  studies  suggest  that  99%  of  dissolved  iron  (i.e.  that  which  passes  through 
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Figure  1-4:  [Fe]  (<  0.4  /xm)  at  the  surface  (top  panel),  at  1000-2000m  (middle  panel) 
and  depths  >  2000m  (bottom  panel).  Data  sources:  Boyle  (unpublished);  de  Baar  et 
al,  1999;  de  Jong  et  al. ,  1998;  Johnson  et  al .,  1997  and  references  therein;  Measures 
and  Vink,  2001;  Nishioka  et  al., 2001;  Obata  et  al,  1993;  Powell  and  Donat,  2001; 
Rue  and  Bruland,  1995;  Sedwick  et  al,  2000;  Sohrin  et  al,  2000;  Takeda  et  al,  1995; 
Wu  and  Luther,  1994;  Wu  et  al,  2001 
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a  0.4  fim  filter)  is  bound  to  organic  ligands  throughout  the  world’s  oceans  (Gledhill 
and  van  den  Berg,  1994;  Rue  and  Bruland,  1995;  van  den  Berg,  1995;  Wu  and  Luther, 
1995;  Rue  and  Bruland,  1997;  Gledhill  et  al.,  1998;  Nolting  et  al.,  1998;  Witter  and 
Luther,  1998;  Witter  et  al.,  2000;  Boye  et  al.,  2001;  Powell  and  Donat,  2001).  The 
reaction  between  Fe'  and  an  Fe-binding  organic  ligand  (L'),  a  molecule  with  potential 
binding  sites,  is: 


Fe'  +  L'^FeL  (1.1) 

where  L'  is  the  Fe-binding  organic  ligand.  The  thermodynamic  equilibrium  can 
be  expressed  as: 


KFeL  =  [FeL)/[Fe'][L' )  (1.2) 

The  source(s),  sink(s)  and  chemical  characterization  of  the  ligand  is  not  well 
known.  Estimates  of  the  concentration  of  ligand  range  between  0.5-6  nM  and  most 
studies  suggest  only  one  class  of  active  organic  ligand,  but  two  studies  (Rue  and 
Bruland,  1997;  Nolting  et  al.,  1998)  have  inferred  two  ligand  classes  in  the  North 
Pacific  and  the  Pacific  sector  of  the  Southern  Ocean.  Vertical  ligand  profiles  appear 
nutrient-like  with  ligand  concentration  below  1000m  remaining  relatively  constant. 

The  estimated  conditional  stability  constant  of  the  ligand(s)  ( KFeL )  ranges  be¬ 
tween  109  8  M-1  and  1014  3  M_1  without  any  clear  regional  pattern.  Conditional  sta¬ 
bility  constants  are  estimated  using  competitive  ligand  exchange  adsorptive  cathodic 
stripping  voltammetry  (CLE-ACSV).  A  known  concentration  of  a  well  characterized 
and  purified  synthetic  ligand  is  added  to  a  series  of  field  samples  containing  natu¬ 
ral  ligands  and  a  range  of  added  metal.  Once  the  synthetic  ligand  has  equilibrated, 
ACSV  is  used  to  measure  the  concentration  of  the  metal  complexed  to  the  added 
ligand  as  a  function  of  total  metal  in  solution. 

Traditionally,  only  the  free  metal  is  considered  to  be  bioavailable  (Anderson  and 
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Figure  1-5:  Schematic  diagram  of  iron  cycling  in  the  global  ocean.  Figure  courtesy 
of  Bridget  Bergquist. 

Morel,  1978).  Since  essentially  the  total  Fe  pool  is  complexed,  the  complexed  Fe  may 
actually  be  bioavailable.  Furthermore,  extracting  iron-binding  compounds  from  sea¬ 
water  and  characterizing  certain  functional  groups  within  the  compound,  Macrellis  et 
al.  (2001)  found  that  functional  groups  known  to  be  present  in  marine  and  terrestrial 
siderophores  are  present  in  the  marine  environment.  This  suggests  that  the  ligands 
are  produced  biologically  by  phytoplankton  to  aid  in  the  uptake  of  Fe  from  seawater. 

The  various  processes  that  iron  is  subjected  to,  which  my  model  aims  to  mecha¬ 
nistically  describe,  are  summarized  in  a  schematic  diagram  (Figure  2-1)  and  will  be 
described  in  more  detail  later. 
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1.2  Prior  Modeling  Studies 


Bruland  et  al.  (1994)  and  Boyle  (1997)  suggested  that  the  variable  aeolian  input  into 
each  ocean  basin  coupled  with  the  biological  processes  of  uptake/regeneration  and  the 
metallic  property  of  scavenging  could  explain  the  profile  of  Fe  in  the  world’s  oceans. 
Johnson  et  al.  (1997)  suggest  that  iron’s  complexation  to  an  organic  ligand  must  be 
controlling  the  deep  water  distribution. 

Lefevre  and  Watson  (1999)  and  Archer  and  Johnson  (2000)  developed  models  to 
examine  possible  controls  on  deep  water  Fe  gradients.  Lefevre  and  Watson  (1999)  use 
a  ten  box  representation  of  the  ocean  adapted  from  the  PANDORA  model  (Broecker 
and  Peng,  1986, 1987),  parameterizing  scavenging  of  iron  onto  particles,  as  well  as  its 
biological  consumption,  remineralization  and  aeolian  deposition.  With  a  scavenging 
rate  of  0.005  yr-1  and  solubility  of  the  aeolian  iron  supply  of  2%,  their  model  was  able 
to  reproduce  the  broad  features  of  the  deep  water  iron  gradients,  although  the  model 
yielded  higher  absolute  concentrations  in  the  deep  Atlantic  (1.6  nM)  than  observed. 
Lefevre  and  Watson  (1999)  also  introduced  a  parameterization  of  complexation  and 
scavenging,  which  assumed  an  effective  solubility  of  iron  (F esoi)  to  represent  the  iron 
complexed  by  a  ligand,  having  a  uniform  oceanic  concentration  of  0.6nM.  Iron  loss 
was  parameterized  as  a  damping  toward  the  effective  solubility,  —k  *  ([Fe]  —  Fesoi), 
with  a  timescale  of  1  /k  =  100  years,  assuming  that  only  iron  which  exceeds  the 
effective  solubility  (i.e.  is  not  bound  to  the  ligand)  can  be  scavenged  from  the  water 
column.  Implicitly,  complexation  to  the  ligand  is  assumed  to  be  very  rapid.  In 
this  model,  deep  water  concentrations  in  the  Antarctic  and  the  Indo-Pacific  were 
approximately  0.6  nM,  and  the  deep  Atlantic  is  somewhat  higher  with  a  concentration 
of  0.92  nM.  The  concentrations  are  higher  than  currently  observed  in  the  Atlantic  and 
Indo-Pacific.  However,  at  the  time  of  that  study,  this  model  seemed  more  consistent 
with  the  available  evidence,  which  was  interpreted  to  show  that  the  deep  water  iron 
concentration  was  uniform  in  all  basins.  There  were  fewer  measurements  and  none 
from  the  Southern  Ocean  at  that  time. 


32 


Archer  and  Johnson  (2000),  using  a  three-dimensional  global  circulation  and  bio¬ 
geochemistry  model  examined  three  parameterizations  of  iron  cycling:  (A)  scavenging 
only,  (B)  complexation  with  one  ligand  and  (C)  complexation  with  two  (strong  and 
weak)  ligands.  In  the  first  case,  using  a  slow  scavenging  rate  (1.6*10~3  yr_1),  the  deep 
water  distribution  reflects  that  of  a  typical  nutrient.  In  the  second  case,  representing 
complexation  with  a  very  strong  ligand  (K=1.2*1013)  of  uniform  concentration  (0.6 
nM)  results  in  a  uniform  deep  water  Fe  distribution,  consistent  with  the  observations 
and  their  interpretation  at  the  time.  In  the  third  case,  Archer  and  Johnson  (2000) 
apply  the  profile  of  two  iron-binding  ligands,  a  strong  ligand  (K=1.2*1013  M-1)  in 
the  upper  500  meters  with  a  maximum  concentration  of  0.5  nM  and  a  weaker  ligand 
(K=3*10h  M-1)  with  concentrations  ranging  between  1.5-2. 5  nM  from  the  surface  to 
depth,  as  measured  by  Rue  and  Bruland  (1995)  in  the  North  Pacific.  This  model  sim¬ 
ulation  results  in  roughly  uniform  deep  water  [Fe]  in  the  Atlantic  and  Pacific  basins. 
In  this  scenario,  excess  iron  at  the  surface  not  utilized  biologically  was  removed  from 
the  system, contributed  to  an  unidentified  process. 

In  the  light  of  new  data,  the  first  goal  of  this  thesis  is  to  adapt  and  constrain 
parameterizations  of  iron  cycling  in  the  deep  ocean.  I  test  three  parameterizations  of 
iron  cycling:  (A)  net  scavenging  onto  particles,  (B)  explicit  representation  of  scav¬ 
enging  and  desorption  to  and  from  particles,  and  (C)  scavenging  and  complexation 
by  a  single  ligand  with  an  imposed  total  concentration.  Here,  I  test  these  parame¬ 
terizations  within  the  context  of  a  computationally  economical,  but  highly  idealized, 
six-box  model  and  a  more  sophisticated  three-dimensional  ocean  general  circulation 
model  (GCM)  to  consider  which  of  these  parameterizations  most  successfully  captures 
the  observed  system  as  present  observations  show. 
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1.3  Increased  Aeolian  Supply 


There  is  particular  interest  in  the  role  of  iron  and  aeolian  dust  supply  in  the  mod¬ 
ulation  of  the  surface  nutrient  concentration  of  the  Southern  Ocean  and  the  ocean’s 
biological  pumps  of  carbon  (Martin,  1990).  Evidence  from  ice  cores  (Petit  et  al,  1999) 
and  suggestions  from  numerical  models  (Mahowald  et  al.,  1999)  indicates  an  increased 
aeolian  supply  of  iron  throughout  the  oceans  during  periods  of  glaciation.  Numeri¬ 
cal  ocean  models  have  been  used  to  explore  the  implications  for  the  carbon  cycle  by 
examining  the  response  to  imposed  surface  nutrient  drawdown  (e.g.  Sarmiento  and 
Orr,  1991).  Watson  et  al.  (2000)  used  a  simplified  ocean  biogeochemistry  model  with 
explicit  representation  of  iron  cycling,  forced  with  glacial-interglacial  cycles  in  South¬ 
ern  Ocean  iron  deposition  derived  from  ice  core  dust  records.  Their  study  suggests 
that  a  significant  fraction  of  the  observed  glacial-interglacial  change  in  atmospheric 
CO2  may  be  accounted  for  in  this  way.  Deep  water  iron  cycling  is  represented  in  that 
model  as  a  particulate  scavenging  process.  In  contrast,  Lefevre  and  Watson  (1999) 
found  it  necessary  to  increase  dust  flux  globally  by  a  factor  of  10  in  order  to  draw¬ 
down  modeled  pC02  by  50  //atm.  Archer  and  Johnson  (2000)  examined  the  response 
of  surface  phosphate  loading  to  increased  aeolian  dust  supply  in  their  global,  three- 
dimensional  model  in  which  deep  ocean  iron  cycling  is  represented  as  a  combination 
of  complexation  to  organic  ligands  and  scavenging  by  particles.  The  authors  suggest 
that  a  significant  drawdown  of  the  surface  macronutrients  might  be  achieved  with 
high  ligand  concentrations. 

The  nature  of  the  parameterization  of  iron,  as  well  as  the  distribution  of  and 
amount  of  dust  flux  increase,  appears  to  affect  the  sensitivity  of  pC02  drawdown  and 
surface  PO4  draw'down  in  the  Southern  Ocean,  but  varying  results  are  drawn  from  a 
diverse  suite  of  models  and  experiments.  The  second  goal  of  this  thesis  is,  therefore,  to 
examine  these  issues  through  a  set  of  sensitivity  experiments  in  a  common  framework 
using  the  idealized  six-box  model.  Further,  I  explore  the  response  of  surface  P04  in 
the  HNLC  regions  to  increased  dust  flux  in  the  three-dimenensional  GCM. 
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1.4  Summary 


It  is  clear  that  iron  plays  an  important  role  in  controlling  the  efficiency  of  the  biological 
pump  and  the  ultimate  drawdown  of  CO2  by  the  ocean.  Yet,  there  are  many  open 
questions  regarding  our  understanding  of  the  marine  iron  cycle  and  how  to  represent 
it  in  ocean  models.  This  thesis  develops  an  iron  biogeochemistry  model  that  draws 
on  previously  published  models  (Archer  and  Johnson,  2000;  Lefevre  and  Watson, 
1999),  but  incorporates  knowledge  from  more  recent  observations  and  experimental 
evidence.  I  incorporate  my  iron  biogeochemical  model  into  a  simple  multi-dimensional 
box  model  and  an  ocean  three-dimensional  general  circulation  model,  to  test  our 
understanding  of  controls  on  deep  water  iron  distribution.  Furthermore,  I  test  the 
effect  that  increased  dust  flux  has  on  surface  P04  drawdown  in  the  Southern  Ocean 
to  place  constraints  on  the  role  iron  may  play  in  glacial-interglacial  cycles. 

In  Chapter  2,  three  parameterizations  of  iron  cycling  are  presented:  a  net  scav¬ 
enging  case,  scavenging-desorption  case  and  complexation-scavenging  case.  Embed¬ 
ded  into  a  simple  six-box  model,  each  parameterization  is  able  to  reproduce  the 
observed  deep  water  gradients.  Results  of  sensitivity  studies  to  various  rate  constants 
and  parameterizations  are  discussed. 

The  sensitivity  of  the  three  iron  parameterizations  to  increased  dust  flux,  as  the 
earth  likely  experienced  during  the  Last  Glacial  Maximum  (LGM),  are  presented  in 
Chapter  3.  While  each  of  the  parameterizations  are  able  to  reproduce,  in  a  broad 
sense,  the  modern  observed  deep  water  [Fe],  their  response  to  increased  dust  flux 
varies  significantly. 

I  add  the  iron  biogeochemical  parameterizations  to  an  ocean  general  circulation 
model  in  Chapter  4.  While  all  three  parameterizations  successfully  reproduced 
observed  deep  water  iron  gradients  in  the  box  model  simulations,  the  ocean  general 
circulation  model,  with  a  more  sophisticated  representation  of  physics,  identifies  some 
differences.  Only  the  complexation  case  is  able  to  easily  capture  the  observed  Fe 
distribution. 
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Since  the  complexation  case  is  the  most  mechanistic  of  the  three  parameterizations 
and  is  best  able  to  reproduce  modern  global  iron  distributions,  I  examine  the  results 
from  this  parameterization  in  greater  depth  in  Chapter  5.  I  compare  modeled  profiles 
to  observations  and  find  good  agreement,  except  at  the  surface  in  high  flux  regions 
such  as  the  North  Atlantic.  I  calculate  the  global  ocean  residence  time  for  Fe  of  ~285 
years,  indicating  that  transport  must  play  a  key  role  in  controlling  Fe  distributions.  I 
also  find  that  globally,  upwelling  accounts  for  40%  of  ’new  iron’  reaching  the  euphotic 
zone.  Lastly,  I  define  a  tracer,  Fe*,  that  indicates  the  degree  to  which  a  water  mass 
is  iron  limited. 

In  Chapter  6,  I  examine  the  response  of  surface  PO4  in  HNLC  regions  to  esti¬ 
mated  LGM  dust  fluxes  (Mahowald  et  al.}  1999)  using  the  ocean  general  circulation 
model.  While  surface  P04  is  depleted  considerably,  there  is  still  excess  PO4  at  the 
surface,  despite  a  2.5-fold  global  increase  in  dust  flux  over  modern  day  estimates. 

I  summarize  the  major  findings  of  the  thesis  in  Chapter  7.  Additionally,  I 
recommend  areas  of  further  study  that  would  be  helpful  from  a  modeling  perspective, 
as  well  as  highlight  additional  research  questions  this  thesis  has  brought  to  light. 
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Chapter  2 


Explorations  of  Biogeochemical 
Iron  Cycling  using  a  Multi-box 
Model 


New  data  from  the  Pacific  and  Southern  Oceans  (de  Baar  et  al ,  1999;  Powell  and 
Donat,  2001;  Sohrin  et  al,  2000;  Wu  et  al,  2001;  Boyle  et  al,  in  prep.)  show  the 
distribution  of  dissolved  iron  in  the  deep  ocean  differs  significantly  from  the  uniformity 
that  previous  models,  such  as  those  of  Lefevre  and  Watson  (1999)  and  Archer  and 
Johnson  (2000)  sought  to  reproduce  and  understand.  Concentrations  now  appear  to 
be  lowest  in  the  Southern  Ocean  and  highest  in  the  Atlantic  basin,  with  a  range  of 
values  in  the  Pacific.  These  models  and  parameterization  need  to  be  revisited  in  the 
light  of  the  new  data. 

In  addition,  recent  measurements  indicate  a  range  in  the  strength  of  the  condi¬ 
tional  stability  constant  and  the  presence  of  a  significant  amount  of  free  ligand.  The 
Lefevre  and  Watson  (1999)  model  does  not  account  for  these  observations.  Archer 
and  Johnson’s  (2000)  model  adds  a  weaker  ligand  in  their  two  ligand  scenario,  but 
still  has  a  strong  ligand  at  the  surface,  which  requires  the  bioavailability  of  Fe  reach¬ 
ing  the  high  latitudes  to  be  reduced  relative  to  the  rest  of  the  model-domain  to  keep 
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surface  [PO4]  high  in  the  high  latitude  surface  waters. 

Here  I  aim  to  build  on  these  previous  studies  and  use  more  recent  data  to  adapt 
and  constrain  the  parameterizations  of  iron  cycling  in  the  deep  oceans.  I  will  also 
explore  the  implications  for  our  understanding  of  the  global  nutrient  and  carbon  cycle. 
In  order  to  allow  significant  exploration  of  parameter  space,  I  use  a  computationally 
economical,  though  highly  idealized,  six-box  model  of  ocean  biogeochemistry  similar 
to  that  of  Lefevre  and  Watson  (1999)  and  Broecker  and  Peng  (1986,  1987).  I  test  three 
parameterizations  of  iron  cycling  with  the  model:  (I)  net  scavenging  onto  particles  (II) 
explicit  representation  of  scavenging  and  desorption  to  and  from  particles,  and  (III) 
scavenging  and  complexation.  I  will  show  that,  provided  appropriate  parameter  values 
are  chosen,  each  of  these  parameterizations  can  reproduce  the  broad  characteristics  of 
the  presently  observed  deep  ocean  dissolved  (<  0.4  //m)  iron  distribution.  However, 
only  the  third  case  concurrently  reproduces  the  observed  deep  water  gradients  and 
the  speciation  of  iron  and  ligand. 

In  the  following  sections  I  will  outline  the  structure  and  mechanics  of  the  six-box 
ocean  biogeochemistry  model  and  discuss  model  results  and  sensitivities  for  each  of 
the  three  iron  parameterizations  outlined  above. 


2.1  Global  Ocean  Biogeochemistry  Model 

I  use  a  six  box  model  (Figure  2-1)  similar  in  construction  to  Broecker  and  Peng’s 
(1986,  1987)  Pandora  Model,  representing  the  surface  and  deep  waters  of  the  At¬ 
lantic,  Antarctic,  and  Indo-Pacific  basins.  Each  basin  is  divided  into  two  layers,  a 
100  meter  surface  layer  where  biological  uptake  of  nutrients  occurs,  and  a  deep  layer. 
Broecker  and  Peng  (1987)  tuned  volume  transports  to  optimize  the  modeled  14C  dis¬ 
tribution.  I  recognize  that  such  highly  idealized  models  have  limitations,  particularly 
for  quantitative  assessments  (Broecker  et  al.,  1999;  Archer  et  al.,  2000;  Follows  et  al., 
2002)  but  they  do  provide  a  useful  framework  in  which  to  develop  clear,  qualitative 
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Figure  2-1:  Schematic  Diagram  of  box  model  adapted  from  Broecker  and  Peng  (1986, 
1987).  The  arrows  represent  volume  transport  (Sv). 

understanding  and  preliminary  sensitivity  studies  of  unconstrained  rates  of  various 
processes. 

2.1.1  Representation  of  Macro- nutrient  Cycling 

The  tracers  explicitly  carried  in  my  model  are  phosphate  (PO4),  dissolved  organic 
phosphorus  (DOP),  total  dissolved  iron  (FeT),  and  particulate  inorganic  iron  (FeP). 
Biological  uptake  and  regeneration  are  indexed  to  phosphorus.  I  illustrate  the  me¬ 
chanics  of  the  model’s  phosphorus  cycle  with  the  prognostic  equations  for  phosphate 
(P04)  and  dissolved  organic  phosphate  (DOP)  for  the  surface  and  deep  Atlantic 
(boxes  i  and  ii  in  Figure  2-1).  For  the  surface: 


(2.1) 

(2.2) 

(2.3) 
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Superscript  numerals  indicate  the  relevant  model  reservoir.  The  first  term  on 
the  right  of  (1)  indicates  transport  by  the  model’s  circulation,  the  second  represents 
biological  export,  and  the  third  term  the  remineralization  of  DOP.  Biological  uptake 
and  export  are  limited  by  light,  phosphate  and  iron  (2.3).  In  conditions  where  Fe  and 
light  are  replete,  I  assume  surface  P04  to  be  the  limiting  nutrient  which  is  exported 
with  a  characteristic  timescale,  1///,  of  about  1  month.  Iron  limitation  is  represented 
by  Michaelis-Menten  kinetics.  The  half  saturation  constant  (Ks)  is  globally  uniform 
but  is  adjusted,  within  the  range  of  measured  values  (Price  et  al.t  1994;  Fitzwater 
et  al.,  1996),  to  optimize  the  modeled  surface  [P04]  and  [Fer]  distributions.  For  the 
deep: 


dPO y 

=  -«  •  VPO?  +  F(1  -  fDOpi)  +  A  DOP”  (2.4) 

dDOPH 

— - —  =  -u  ■  xjDOPn  -  A DOPn  (2.5) 

In  (2.4),  the  first  term  on  the  right  represents  transport,  the  second  the  reminer¬ 
alization  of  sinking  particulate  matter,  and  the  third  the  remineralization  of  DOP. 

Two-thirds  of  the  exported  nutrient  (fDOp)  enters  the  surface  dissolved  organic 
phosphorus  (DOP)  pool,  while  one-third  is  rapidly  exported  as  particulate  P  to  the 
deep  P04  pool  (Yamanaka  and  Tajika,  1997).  The  imposed  timescale  for  remineral¬ 
ization  of  DOP,  (1/A),  is  6  months. 

2.1.2  Iron  Cycling 

The  aeolian  source  of  iron  is  prescribed,  while  the  loss  of  iron  due  to  scavenging 
and  iron’s  role  in  the  biological  cycle  are  modeled  explicitly.  Total  dissolved  iron 
(Fer),  and  particulate  inorganic  iron  (FeP)  are  prognostic  tracers  of  the  model.  The 
following  equations  below  describe  the  iron  cycle  for  the  surface  and  deep  Atlantic 
(boxes  i  and  ii  in  Figure  2-1).  The  equations  for  the  other  basins  are  similar.  For  the 
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surface: 


aF\n  -  u  •  V^4  ~  ri?Fe  +  XDOPtRpe  +  4e  (2-6) 

-4e  -  WS§~ZF^P  (2’7) 

The  first  term  on  the  right  in  (2.6)  represents  the  aeolian  source,  the  second  term 
ocean  transport  of  total  iron,  the  third  term  biological  utilization.  Remineralization 
of  DOM  is  represented  by  the  fourth  term  on  the  right  multiplied  by  the  Redfield 
ratio  (RFe)  between  P04  and  Fe.  The  fifth  term,  4e,  represents  the  interactions 
with  particles  or  ligands,  which  differ  between  each  of  the  three  parameterizations 
and  will  be  described  in  more  detail  below. 

Aeolian  deposition  (F,„)  is  the  source  of  iron  to  the  model  ocean.  Iron  deposition 
data  from  Gao  et  a/.(2001),  Duce  and  Tindale  (1991)  and  Jickells  and  Spokes  (2001) 
were  used  to  estimate  the  source  to  the  surface  waters  of  each  basin.  Table  2.1 
summarizes  the  various  datasets  and  the  values  used.  The  solubility  of  Fe  aerosols 
(a)  in  seawater  is  not  well  known,  although  recent  studies  suggest  it  may  be  below  5% 
(Jickells  and  Spokes,  2001;  Spokes  and  Jickells,  1996).  Based  on  results  of  parameter 
space  exploration,  we  use  a  —  0.01  for  the  models  discussed  here. 

Iron  is  biologically  utilized  in  proportion  to  P04  with  a  fixed  Fe:C  ratio  (R Fe)  and 
a  C:P  Redfield  ratio  of  16:1.  Sunda.  and  Huntsman  (1995)  have  published  estimates 
for  the  Fe:C  ratio  that  indicate  marine  phytoplankton  decrease  their  cellular  iron  re¬ 
quirement  to  optimize  growth  in  Fe-stressed  environments,  but  I  have  not  represented 
this  variability  here,  as  a  clear  relationship  has  not  been  established.  The  Fe:C  ratio 
is  globally  adjusted  to  optimize  surface  [PO4]  and  [Fex]  • 

Evidence  from  Th  isotopes  indicates  that  the  mean  sinking  rate  of  fine  particles  is 
between  500  and  1,000  myr-1  (Cochran  et  al.,  1993).  In  order  to  very  crudely  account 
for  the  different  sinking  rates  of  large  and  small  particles,  we  have  assumed  that  10% 


dFe? 

dt 

dFelP 

dt 
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Table  2.1:  Aeolian  Fe  dust  data  (g  Fe  yr-1). 

Basin _ Gao  Duce/Tindale  Jickells/Spokes  Model  values 


Atlantic 

7.73 

8.54 

6.46 

6.46 

Southern  Ocean 

0.071 

- 

- 

0.071 

Indo-Pacific 

5.71 

23.5 

10.17 

10.17 

Data  sets:  Duce 

and  Tindale,  1991; 

Gao  et  al.,  2001; 

and  Jickells  and  Spokes,  2001 

of  particles  are  large  with  a  sinking  rate  of  20,000  myr-1  and  90%  are  small  particles 
with  a  sinking  rate  of  1,000  myr-1,  yielding  an  average  sinking  rate  (Ws)  of  2,900 
myr-1. 

The  deep  equations  for  iron  are: 


dF 4 
dt 


=  -u  •  \jF<%  +  ri?Fe(l  -  fD0P)  +  \DOPHRFe 


+  jf 


Fe 


dFelp 

dt 


-jr 


Fe 


Ws—Fep 

oz 


(2.8) 

(2.9) 


I  examine  three  different  parameterizations  for  the  geochemical  processes:  (I)  net 
scavenging  (II)  scavenging  and  desorption,  and  (III)  scavenging  and  complexation.  In 
case  I  and  II,  I  do  not  differentiate  between  complexcd  and  free  iron  and  assume  that 
the  total  iron  pool  is  subjected  to  all  geochemical  processes.  In  case  III,  I  explicitly 
model  complexation  and  differentiate  between  free  iron  and  complexed  iron. 


2.2  Model  Results 

While  I  will  focus  on  the  iron  distribution  in  this  discussion,  the  phosphate  distri¬ 
bution,  which  is  explicitly  controlled  by  iron  limitation  here,  also  provides  a  consis¬ 
tency  check  on  the  model.  For  solutions  in  which  iron  distribution  is  reasonable,  the 
phosphate  distributions  are  in  good  agreement  with  observations.  Surface  [PO4]  is 
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Symbol  Definition 


Table  2.2:  Model  Parameters 

Value 


R Fe  Fe:P  ratio 

f nop  fraction  of  DOP 

h  Depth  of  surface  box 

H  Depth  of  deep  box 

Fin  Aeolian  Deposition  rate 

a  Fe  dust  solubility 

/j,  Biological  uptake  rate 

A  Remineralization  rate 

k.sc  Scavenging  rate 

kb  Backscavenging  rate 

Ws  Particle  sinking  velocity 

Ks  Iron  half  saturation  constant 

Kpei  Ligand  conditional  stability  constant 


adjustable  parameter 
0.67 
100m 
3900m 
see  Table  1 
1% 

1  month-1 

0.5  yr— 1 

variable 

variable 

2,900  myr-1 

adjustable  parameter 

variable 


elevated  in  the  Southern  Ocean  box,  depleted  in  the  Atlantic  box  and  intermediate 
in  the  Indo-Pacific  box.  Deep  [PO4]  increases  from  the  Atlantic  to  the  Indo-Pacific. 


2.2.1  Case  I:  Net  Scavenging  Model 

Boyle  (1997)  suggested  that  the  deep  water  distribution  of  iron  may  be  modeled  using 
simple  parameterizations  of  aeolian  deposition,  biological  uptake  and  remineralization 
of  organic  matter,  and  a  representation  of  net  scavenging  to  particles.  Such  a  model 
is  highly  idealized,  and  does  not  attempt  to  explicitly  represent  the  details  of  the 
biogeochemical  processes,  but  it  could  be  the  simplest  viable  prognostic  model  for 
iron  in  the  ocean.  It  has  only  one  adjustable  parameter  and  does  not  attempt  to 
describe  poorly  understood  details  of  the  biogeochemical  processes. 

Here,  I  examine  whether  this  parameterization  can  reproduce  the  broad  basin  to 
basin  and  surface  to  deep  ocean  observed  gradients  of  dissolved  iron.  In  this  formula¬ 
tion  FeT  is  scavenged  and  utilized  biologically.  This  parameterization  is  conceptually 
similar  to  the  no-ligand  model  of  Lefevre  and  Watson  (1999).  I  impose  the  regional 
variations  in  aeolian  supply  and  examine  the  sensitivity  of  the  dissolved  iron  distri¬ 
bution  to  the  net  scavenging  rate.  In  this  case  the  loss  of  iron  is  modeled  simply 
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Figure  2-2:  Sensitivity  of  deep  [Fer]  to  scavenging  rates.  For  slow  scavenging  rates, 
(knetsc  <  001),  the  FeT  distribution  is  nutrient-like.  For  intermediate  scavenging  rates, 
.004<  knetsc  <-006,  the  observed  gradients  are  reproduced.  For  kneisc>.006,  the  sense 
of  gradient  is  reproduced,  although  mean  concentrations  are  lower  than  observed. 

as  a  first-order  scavenging  process,  limited  by  the  dissolved  free  Fe  concentration. 
Scavenged  iron  is  transfered  to  the  particulate  pool,  Fep,  with  rate  constant  —knetSc, 
and  is  stripped  from  the  water  column  as  the  particles  sink.  For  the  net  scavenging 
(and  scavenging-desorption)  case  Fe:C  is  25  jumobl  mol.  Here,  then 

J Fe  —  knc,tScFeT-  (2.10) 

Figure  2-2  shows  the  deep  ocean,  dissolved  iron  concentrations  in  each  of  the 
three  modeled  regions  (Atlantic  basin,  Southern  Ocean,  Pacific  basin)  as  a  function 
of  the  net  scavenging  rate.  Each  cluster  of  three  bars  represents  the  solution  of  the 
model  at  a  particular  value  of  scavenging  rate.  The  relative  lengths  of  the  three  bars 
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reflect  the  basin  to  basin  gradients  of  deep  iron  in  each  solution.  In  the  case  of  a 
slow  net-scavenging  rate  (knetsc  =  0.001  yr-1),  the  deep  water  distribution  is  that  of  a 
typical  nutrient  with  the  deep  Indo-Pacific  iron  concentration  greater  than  the  deep 
Southern  Ocean  which  is  greater  than  the  deep  Atlantic.  The  result  is  unsurprising, 
but  the  gradients  are  not  as  observed.  For  stronger  scavenging,  knetsc  >-004  yr-1, 
the  observed  deep  water  Fe  gradients  (Atl>Indo-Pacific>Southern  Ocean)  are  repro¬ 
duced.  However,  when  knetsc  >  006  yr-1,  though  the  inter-basin  gradients  remain  of 
the  correct  sign,  the  mean  ocean  deep  water  [Fe^]  becomes  much  too  low. 

This  simple  model,  representing  the  basin  variations  of  the  aeolian  supply  and 
a  uniform,  net  scavenging  rate  can  reproduce  the  unique  deep  water  iron  signature 
provided  that  0.004  yr-1<  k$c  <  0.006  yr-1.  This  is  consistent  with  the  previous 
study  of  Lefevre  and  Watson  (1999). 

2.2.2  Case  II:  Scavenging-Desorption  Model 

While  the  highly  simplified  model  of  Case  I  can  reproduce  the  broad,  basin-to-basin 
gradients  of  the  dissolved  iron  distribution,  it  does  not  resolve  the  biogeochemical 
processes  at  work.  In  Cases  II  and  III,  I  introduce  more  detailed  parameterizations 
which  attempt  to  represent  processes  known  to  be,  or  likely  to  be,  at  work  in  the 
ocean.  I  ask  if  these  more  detailed  models  can  reproduce  the  observations  and,  if  so, 
what  constraints  can  be  placed  on  system  parameters  by  the  observations? 

Thorium  is  produced  in  the  ocean  by  radio-decay  and  is  subsequently  scavenged 
out  of  the  water  column  by  sinking  particles.  Bacon  and  Anderson  (1982),  using 
oceanic  observations  of  thorium  isotopes,  have  estimated  a  scavenging  rate  between 
0.2-1.28  yr-1.  This  is  much  faster  than  the  net  scavenging  rate  for  iron  implied  in 
our  model  (case  I),  but  does  not  represent  the  net  scavenging  rate  for  Th.  Bacon 
and  Anderson  (1982)  suggest  that  scavenged  Th  is  also  desorbed  from  particles,  i.e. 
released  back  to  the  water  column,  and  also  infer  from  data  a  rate  at  which  this 
occurs.  Since  Fe  and  Th  have  similar  metallic  properties,  I  consider  it  likely  that  iron 
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Figure  2-3:  Schematic  description  of  the  scavenging  and  desorption  model.  Desorption 
is  treated  as  a  first  order  process  dependent  on  the  particulate  iron  concentration  and 
transfers  particulate  Fe  to  the  dissolved  pool.  Scavenging  is  modeled  as  a  first  order 
process  dependent  on  the  dissolved  Fe  concentration.  Scavenged  iron  can  be  lost  from 
the  ocean,  ultimately  balancing  the  aeolian  sink. 


may  experience  a  similar  dynamic  interplay  of  scavenging  and  desorption  to  and  from 
particles. 

To  address  this  possibility  in  case  II,  I  parameterize  the  interactions  of  iron  with 
particles  in  the  deep  water  as  a  cycle  of  rapid  scavenging  and  desorption,  which  may 
result  in  a  slow,  net  scavenging  consistent  with  the  observed  distribution  and  Case  I 
above  (Figure  2-3).  In  this  case 


JFe  —  —kscFe.T  +  kiFep  (2-11) 

Here  —ksc  is  the  scavenging  rate.  Scavenging  is  proportional  to  the  availability  of 
dissolved  iron;  kb  is  the  desorption  rate,  and  desorption  is  proportional  to  partic¬ 
ulate  iron.  Figure  2-4  shows  the  deep  water,  dissolved  iron  concentration  in  each 
of  the  model  regions  as  a  function  of  scavenging  rates  ranging  between  0.1-1  yr-1 
and  desorption  rates  between  20-100  yr-1.  When  the  ratio  of  desorption/scavenging 
~150-170,  this  model  is  able  to  broadly  reproduce  the  observed  global  deep  water  Fe 
gradients  and  concentrations  (dashed  contours). 
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Figure  2-4:  Scavenging-desorption  model:  Sensitivity  of  deep  [Fer]  (nM)  to 
scavenging-desorption  rate  constants.  [Fey]  as  a  function  of  scavenging  (ksc,  yr-1) 
and  desorption  (kb,  yr-1)  for  A)  Atlantic,  B)  Southern  Ocean  and  C)  Indo-Pacific 
basin.  The  dashed  contours  indicate  the  average  observed  [Fey]  for  each  basin.  The 
optimal  solution  is  for  kb/ksc  ~150-170  . 
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For  thorium,  the  desorption  to  scavenging  ratio  is  calculated  to  be  an  order  of 
magnitude  smaller.  I  might  interpret  these  model  results  to  suggest  that  iron  and 
thorium  may  behave  in  a  similar  manner,  but  have  different  desorption  to  scavenging 
ratios,  implying  that  Th  is  much  more  particle  reactive  than  Fe.  On  the  other  hand, 
there  are  other  processes  which  may  be  significant  for  iron  and  which  I  should  include 
in  the  model. 

2.2.3  Case  III:  Complexation 

Case  II  again  found  a  plausible  solution  of  the  model  by  representing  iron  as  an  ana¬ 
logue  of  thorium,  provided  appropriate  scavenging  and  desorption  rates  are  applied. 
New  methods  and  observations  of  iron  in  the  ocean  would  be  required  to  directly 
confirm  such  a  mechanism  at  work.  However,  there  is  a  great  deal  of  evidence  that 
another  biogeochemical  process  -  complexation  with  organic  ligands  -  plays  a  signifi¬ 
cant  role  in  the  control  of  deep  water  iron  distributions. 

Observational  evidence  (Gledhill  and  van  den  Berg,  1994;  Rue  and  Bruland,  1995; 
van  den  Berg,  1995;  Wu  and  Luther,  1995;  Rue  and  Bruland,  1997;  Gledhill  et  al., 
1998;  Nolting  et  al.,  1998;  Witter  and  Luther,  1998;  Witter  et  al.,  2000;  Boye  et  al., 
2001;  Powell  and  Donat,  2001)  indicates  that  over  99%  of  “dissolved”  iron  is  bound 
to  a  ligand.  In  this  third  case  I  add  a  mechanistic  description  of  Fe  complexation 
to  our  box  model  (Figure  2-5).  Representations  of  the  effect  of  complexation  have 
been  introduced  in  two  previous  models  (see  Chapter  1).  The  model  applied  here 
is  closely  related  to  the  (second)  model  of  Archer  and  Johnson  (2000),  representing 
complexation  with  a  single  ligand  imposing  [Lx].  In  the  light  of  new  data,  I  explore 
the  sensitivity  of  modeled  deep  water  [Fex]  to  a  wide  range  of  parameter  values.  For 
example,  in  the  Archer  and  Johnson  (2000)  model,  LT  is  fixed  at  0.6  nM,  while  I  test 
the  sensitivity  of  deep  water  [Fer]  to  the  value  of  [Lx],  which  varies  between  0.5-6  nM 
in  the  global  ocean  (see  Chapter  1.1).  I  also  assess  the  sensitivity  of  deep  water  [Fex] 
to  a  range  of  scavenging  rates  and  ligand  strengths,  since  they  arc  poorly  constrained. 
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Figure  2-5:  Schematic  diagram  of  the  complexation  model.  Dissolved  Fe  can  undergo 
two  transformations:  it  can  be  scavenged  or  it  can  be  complexed.  The  box  represents 
the  reaction  Fe'+L'=FeL.  We  assume  that  chemical  forms  within  the  box  (Fe'  and 
FeL)  can  be  utilized  biologically,  but  only  Fe'  can  be  scavenged. 

Here,  dissolved  iron  is  assumed  to  be  the  sum  of  “free”  and  “complexed”  forms: 

FeT  =  Fe!  +  FeL.  (2.12) 

where  FeL  represents  the  iron  complexed  with  an  organic  ligand.  Only  the  free  form 
is  available  for  scavenging  and  hence 

JFe  =  —kscFe'  .  (2.13) 

Since  complexation  occurs  on  the  timescales  of  minutes  to  hours  (Witter  et  al., 2000), 
it  is  assumed  that  the  reaction  goes  to  equilibrium.  I  specify  the  total  ligand  concen¬ 
tration,  L7'=[FeL]+[L'],  and  use  the  equilibrium  relationship  K^/=k///crf=[FeL]/[Fe']  [L/] 
to  determine  the  speciation  of  the  iron.  FeT  is  a  conserved  property  and  is  advected. 
Desorption  from  particles  is  neglected  in  this  case  since  its  impact  is  overwhelmed 
by  the  strong  complexation  reaction.  In  the  complexation-scavenging  case  the  Fe:C 
ratio  is  set  to  15  /rmol:l  mol. 
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Setting  Lx  to  1  nM  in  Figure  2-6, 1  plot  the  relationship  of  the  deep  water  dissolved 
iron  concentration  in  each  basin  to  scavenging  rate,  ranging  between  0.2-1.8  yr-1  and 
conditional  stability  constant,  KFeL,  between  1010  M"1  and  1013  M"1,  reflecting  the 
range  of  values  inferred  from  ocean  observations  (Gledhill  and  van  den  Berg,  1994; 
Rue  and  Bruland,  1995;  van  den  Berg,  1995;  Wu  and  Luther,  1995;  Rue  and  Bruland, 
1997;  Gledhill  et  al,  1998;  Nolting  et  al.,  1998;  Witter  and  Luther,  1998;  Witter  et  al., 
2000;  Boye  et  al.,  2001;  Powell  and  Donat,  2001).  Deep  iron  concentrations  generally 
increase  with  increasing  stability  constant  and  decreasing  scavenging  rate.  Since  only 
the  uncomplexed  form  of  iron  can  be  scavenged,  at  high  scavenging  rates  a  strong 
ligand  is  required  to  maintain  deep  water  ''dissolved”  [Fex]  concentrations  at  observed 
levels,  sequestering  it  in  a  form  which  I  assume  is  not  available  for  scavenging.  At  very 
low  scavenging  rates,  the  sensitivity  to  the  conditional  stability  constant  decreases, 
since  it  is  no  longer  necessary  for  iron  to  be  in  complexed  form  to  remain  in  the 
water  column  for  a  significant  period.  The  sensitivity  to  the  scavenging  constant  is 
weak  when  scavenging  is  strong  because  there  is  very  little  scavengable  iron  and  the 
limiting  process  is  complexation. 

Observations  indicate  that,  while  most  “dissolved”  iron  is  in  complexed  form, 
a  significant  fraction  of  ligand  is  free  (Gledhill  and  van  den  Berg,  1994;  Rue  and 
Bruland,  1995;  van  den  Berg,  1995;  Wu  and  Luther,  1995;  Rue  and  Bruland,  1997; 
Gledhill  et  al,  1998;  Nolting  et  al.,  1998;  Witter  and  Luther,  1998;  Witter  et  al,  2000; 
Boye  et  al,  2001;  Powell  and  Donat,  2001).  This  is  in  contrast  to  the  models  of  Archer 
and  Johnson  (2000)  and  Lefevre  and  Watson  (1999)  where,  due  to  the  low  total  ligand 
concentration  and  high  conditional  stability  constant,  the  dissolved  iron  concentration 
was  about  the  same  as  the  total  ligand  concentration  (0.6  nl VI)  over  much  of  the  ocean. 
This  case,  where  the  ligand  is  saturated,  represents  a  limit  case  of  the  scheme  used 
here.  By  relaxing  these  constraints,  it  is  possible  to  find  a  solution  consistent  with 
the  observed  iron  distribution  which  also  predicts  a  significant  presence  of  free  ligand, 
L'.  Figure  2-7  shows  the  dependency  of  [L']  on  Kpe>L  and  scavenging  rate  constant 
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Figure  2-7:  Complexation  model:  Sensitivity  of  the  free  ligand  concentration  ([I/]) 
(nM)  to  scavenging  rate  (ksc,  yr_1)  and  conditional  stability  constant  (log  Kjv£,)  for 
the  A)  Atlantic,  B)  Southern  Ocean  and  C)  Indo-Pacific  basin  with  [LT]=  1  nM.  As 
log  Kpe>L  increases.  [L]  decreases  due  to  forward  reaction  L'+Fe'=FeL.  As  scavenging 
increases,  [L]  increases,  as  forward  reaction  is  limited  by  Fe,  resulting  in  excess  L. 
Observations  indicate  L y  ranges  between  0.5-6  nM. 
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for  this  model  with  specified  total  ligand  concentration  of  1  nM.  As  the  scavenging 
rate  increases,  the  loss  of  Fe  limits  the  complexation  reaction,  resulting  in  excess  free 
ligand,  [L'].  Comparing  Figure  2-6  and  Figure  2-7,  FeT  and  L  are  inversely  related. 
For  strong  KpeL ,  Fer^FeL,  which  is  the  limit  modeled  by  Archer  and  Johnson  (2000) 
and  implicitly  by  Lefevre  and  Watson  (1999)  . 

Observations  also  indicate  a  significant  variation  in  ligand  concentration  around 
the  ocean  but,  as  yet,  without  a  clearly  emerging  large  scale  pattern  (Gledhill  and 
van  den  Berg,  1994;  Rue  and  Bruland,  1995;  van  den  Berg,  1995;  Wu  and  Luther, 
1995;  Rue  and  Bruland,  1997;  Gledhill  et  al. ,  1998;  Nolting  et  al. ,  1998;  Witter  and 
Luther,  1998;  Witter  et  al.,  2000;  Boye  et  al.,  2001;  Powell  and  Donat,  2001).  Still 
without  introducing  any  spatial  variations  in  the  ligand  concentration,  I  also  illustrate 
the  sensitivity  of  dissolved  iron  and  free  ligand  concentrations  to  the  concentration  of 
total  ligand.  Figure  2-8  and  Figure  2-9  show  the  deep  ocean  iron  concentration  and 
free  ligand  concentration  respectively,  (as  Figures  2-6  and  2-7),  but  with  increased 
total  ligand  concentration,  L t—  4  nM.  For  identical  choices  of  ksc  and  KpeL  with 
increased  total  ligand,  I  find  increased  [Fer].  Hence,  to  the  fit  the  modern  observed 
distribution  with  Ly-=  4  nM,  I  must  adjust  ksc  by  a  factor  of  ~  15-25  times.  However, 
the  sensitivity  pattern  is  the  same  as  the  case  where  L t  is  fixed  to  1  nM. 

The  model  predicts  an  excess  [L']  ranging  from  0.5-3  nM  for  scavenging  rates 
between  0.2-1. 8  yr_1  and  ligand  strengths  ranging  from  log (KFeL)  of  10-13  (Figure 
2-9).  It  suggests  highest  excess  [L']  for  the  Atlantic  basin,  in  broad  agreement  with 
observations  (Gledhill  and  van  den  Berg,  1994;  Rue  and  Bruland,  1995;  van  den  Berg, 
1995;  Wu  and  Luther,  1995;  Rue  and  Bruland,  1997;  Gledhill  et  al,  1998;  Nolting  et 
al.,  1998;  Witter  and  Luther,  1998;  Witter  et  al,  2000;  Boye  et  al,  2001;  and  Powell 
and  Donat,  2001). 


53 


Figure  2-8:  Complexation  model:  Sensitivity  of  modeled  deep  water  [Fe^]  to  scaveng¬ 
ing  rate  ( ksc ,  yr_1)  and  conditional  stability  constant  (log  KFe>L )  for  the  A)  Atlantic, 
B)  Southern  Ocean,  and  C)Indo-Pacific  basin  with  [L^]  =  4  nM.  The  dashed  contour 
represents  the  observed  average  deep  water  [Fer]  for  each  basin. 

2.3  Discussion 

I  have  examined  three  parameterization  of  water  column  iron  biogeochemistry  in 
the  framework  of  an  idealized,  six-box  ocean  biogeochemistry  model.  In  the  light  of 
the  latest  available  observations  of  the  deep  ocean  distribution  of  iron,  an  extremely 
simple  model  which  parameterizes  deep  ocean  biogeochemical  cycling  of  iron  as  a 
first  order  net  scavenging,  is  able  to  capture  the  broad  basin  to  basin  [Fer]  gradients. 
However,  this  parameterization  does  not  explicitly  represent  the  processes  believed 
to  control  the  system.  A  second  parameterization  treated  iron  as  an  analogue  of 
thorium,  with  rapid  scavenging  and  desorption  of  iron  to  and  from  particles.  For  a 
scavenging-desorption  rate  constant  of  ~150,  this  model  can  also  reproduce  the  broad 
features  of  the  large  scale  distribution  of  dissolved  iron. 

In  a  third  parameterization,  following  Archer  and  Johnson  (2000),  I  introduce 
complexation  to  an  organic  ligand.  Sensitivity  studies  showed  that  this  model  can 
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Figure  2-9:  Complexation  model:  Same  as  Figure  2-7,  but  [LT]=  4  nM. 


reproduce  the  large  scale  iron  distribution  for  a  range  of  ligand  strengths,  KpeL  and 
scavenging  rates,  ksc,  (in  appropriate  combination)  and  for  a  range  of  total  ligand 
concentrations,  LT.  The  ligand  parameterization  of  Lefevre  and  Watson  (1999)  and 
Archer  and  Johnson’s  (2000)  complexation  with  one  ligand  case,  with  a  very  strong 
ligand  and  low  total  ligand  concentration,  both  led  to  quite  uniform  deep  ocean  iron 
distributions  and  saturated  ligand.  This  is  a  limit  case  of  the  more  general  model 
presented  here.  The  model  and  recent  observational  data  suggest  that  the  parameter 
choices  of  Archer  and  Johnson’s  (2000)  two  ligand  model,  with  a  very  strong  ligand 
in  the  upper  ocean  resulting  in  fairly  uniform  deep  water  [Fer],  is  at  odds  with  recent 
observational  evidence.  It  would  also  lead  to  high  iron  and  low  phosphorus  concen¬ 
trations  at  the  surface.  To  prevent  the  accumulation  of  iron  in  surface  waters,  Archer 
and  Johnson  (2000)  remove  any  surface  iron  from  the  system  that  is  not  utilized 
biologically,  but  the  process  this  should  represent  is  not  clearly  identified.  Based  on 
the  sensitivity  studies  performed  here  and  recent  observational  data,  I  suggest  that 
a  parameter  regime  with  a  weaker  ligand  and  greater  concentration  of  total  ligand 
may  be  more  realistic.  In  the  latter  case,  the  model  can  reproduce  both  the  deep  iron 
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distribution,  but  also  the  observed  presence  of  significant  amounts  of  free  ligand. 
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Chapter  3 


Sensitivity  of  Surface  Phosphate  to 
Aeolian  Iron  Source  and  Upwelling 
Strength 


A  strong  motivation  for  modeling  the  iron  cycle  is  to  be  able  to  explicitly  describe 
and  explore  the  role  of  iron  in  setting  current,  past  and  future  ocean  distributions  of 
carbon  and  macronutrients.  Of  particular  interest  is  the  possible  impact  and  feed¬ 
backs  of  climate  change  and  the  aeolian  supply  of  iron  to  the  efficiency  of  the  carbon 
pump  in  the  remote  Southern  Ocean.  Martin  (1990)  suggested  increased  dust  flux 
during  the  Last  Glacial  Maximum  (LGM)  could  have  increased  export  production 
in  the  Southern  Ocean  and  decreased  atmospheric  pC02.  While  data  from  ice  cores 
(Petit  et  al,  1999)  and  atmospheric  dust  models  (Mahowald  et  al,  1999)  suggest 
that  aeolian  iron  supply  increased  up  to  twenty  times  in  the  Southern  Ocean  and 
globally  2-5  times  relative  to  present  day,  paleo  productivity  proxies  do  not  suggest 
that  export  production  was  higher  during  the  Last  Glacial  Maximum  (LGM)  in  the 
Southern  Ocean  (Francois  et  al,  1997;  Kumar  et  al,  1995).  Rather  <515N  data  sug¬ 
gests  increased  efficiency  of  nutrient  utilization  in  the  high  latitudes,  perhaps  due  to 
increased  stratification  between  surface  and  deep  waters  leading  to  weaker  vertical 
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exchange  (Francois  et  al.,  1997).  Keeling  and  Visbeck  (2001)  argue  that  the  response 
of  ocean  eddies  to  increased  stratification  would  increase  the  upwelling  rate  of  deep 
waters  in  the  Southern  Ocean.  They  propose  instead  that  enhanced  sea-ice  coverage 
in  the  Southern  Ocean  could  have  prevented  the  outgassing  of  C02. 

Watson  et  al.  (2000)  used  a  simplified  ocean  biogeochemistry  model  with  ex¬ 
plicit  representation  of  iron  cycling,  forced  with  glacial-interglacial  cycles  in  Southern 
Ocean  iron  deposition  derived  from  ice  core  dust  records.  Representing  deep  water 
iron  cycling  as  a  particulate  scavenging  process,  their  study  suggests  that  a  signifi¬ 
cant  fraction  of  the  observed  glacial-interglacial  change  in  atmospheric  C02  may  be 
accounted  for  by  increased  export  production  due  to  the  increased  dust  supply.  This 
is  in  contrast  to  the  results  of  a  study  using  a  six-box  model  in  which  complexa- 
tion  and  scavenging  of  iron  are  parameterized  (Lefevre  and  Watson,  1999).  There, 
despite  increasing  Southern  Ocean  dust  flux  by  a  factor  of  twenty,  pC02  drawdown 
was  by  only  1  //atm.  However,  the  latter  model  did  suggest  a  pC02  drawdown  of 
50//atm  when  global  dust  flux  increased  ten  times.  Lefevre  and  Watson  (1999)  also 
found  very  little  sensitivity  to  the  strength  of  upwelling.  Archer  and  Johnson  (2000) 
also  examined  the  response  of  surface  phosphate  loading  to  increased  aeolian  dust 
supply  in  their  global,  three-dimensional  model  in  which  deep  ocean  iron  cycling  is 
represented  as  a  combination  of  complexation  to  organic  ligands  and  scavenging  by 
particles.  They  show  that  a  significant  drawdown  of  the  surface  macro-nutrients  can 
be  achieved  with  high  ligand  concentrations.  The  nature  of  the  parameterization  of 
iron,  as  well  as  the  distribution  of  and  amount  of  dust  flux  increase,  appears  to  affect 
the  sensitivity  of  pC02  drawdown  and  surface  P04  drawdown  in  the  Southern  Ocean, 
but  varying  results  are  drawn  from  a  diverse  suite  of  models  and  experiments. 
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3.1  Results:  Increased  Aeolian  Flux  Simulations 


Here  I  examine  these  issues  and  apparent  contradictions  through  a  set  of  sensitiv¬ 
ity  experiments  in  a  unified  model  framework,  the  six-box  model  described  in  the 
previous  chapter.  I  present  results  that  illustrate  the  sensitivity  of  surface  phosphate 
drawdown  to  1)  the  magnitude  and  distribution  of  the  dust  supply  and  2)  the  strength 
of  the  vertical  exchange  between  the  Southern  Ocean  surface  and  deep  waters  for  the 
three  parameterizations  (see  Chapter  2)  of  iron.  I  find  considerable  differences  in  the 
response  due  to  the  parameterization  of  iron  and  the  distribution  and  magnitude  of 
the  increased  dust  flux. 

3.1.1  Increased  Dust  Flux:  Southern  Ocean 

I  test  the  sensitivity  of  Southern  Ocean  surface  PO4  to  increased  dust  flux  in  only  the 
Southern  Ocean  and  vertical  exchange  for  each  of  the  parameterizations  as  described 
in  Chapter  2.  In  Figure  3-1,  I  plot  the  Southern  Ocean  surface  [PO4]  from  the  model 
as  a  function  of  increased  aeolian  dust  flux  in  the  Southern  Ocean  only,  relative  to 
today’s,  and  for  several  rates  of  Southern  Ocean  vertical  mass  exchange.  Even  with 
vertical  exchange  weakened  by  50  %  and  and  Southern  Ocean  aeolian  flux  increased 
20  times  relative  to  modern  day  flux  in  the  Southern  Ocean,  surface  PO4  drawdown 
is  only  0.3  pM  for  the  scavenging  and  scavenging-desorption  case  and  0.45  /xM  for  the 
complexation  case.  Yet  the  model  of  Watson  et  al.  (2000)  had  predicted  a  maximum 
drawdown  of  PO4  of  0.6  /xmol  with  a  25  fold  increase  in  the  dust  flux  reaching  the 
Southern  Ocean  above  modern  day  flux.  I  must  increase  Southern  Ocean  dust  flux 
more  than  50  fold  to  observe  a  similar  drawdown  of  PO4  with  my  model  for  all 
three  parameterizations.  While  my  model  results  agree  qualitatively  with  Watson  et 
al.  (2000),  quantitative  differences  in  the  sensitivity  of  surface  PO4  to  increases  in 
dust  flux  can  potentially  be  attributed  to  differences  in  the  parameterization  of  iron. 
Watson  et  al.  (2000)  parameterize  the  loss  of  iron  due  to  particulate  scavenging, 
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Southern  Ocean  Dust  Increase 

Figure  3-1:  Steady  state  surface  Southern  Ocean  P04  sensitivity  to  Southern  Ocean 
only  dust  increase  (absolute  factor)  and  S.  Ocean  overturning  (Sv)  for  A)  net  scav¬ 
enging  case,  B)  scavenging-desorption  case,  and  C)  complexation  case. 


with  the  scavenging  rate  varying  as  a  function  of  particulate  organic  matter  flux.  In 
addition,  the  Fe:C  ratio  is  variable  in  the  Watson  et  al.  (2000)  model.  It  is  interesting 
to  note  that  a  fifty  fold  increase  in  Southern  Ocean  aeolian  flux  is  approximately  equal 
to  estimates  of  current  day  average  aeolian  flux  to  the  Atlantic  basin  per  unit  area. 
My  results  suggest  that  deep  water  transport  of  iron  from  the  other  basins  may  be 
necessary  for  Fe  to  affect  the  drawdown  of  surface  P04. 


3.1.2  Increased  Dust  Flux:  Globally 

Since  increasing  dust  flux  in  only  the  Southern  Ocean  does  not  result  in  significant 
increased  utilization  of  surface  P04,  I  test  the  sensitivity  of  surface  P04  to  a  global 
increase  in  dust  flux.  In  Figure  3-2,  for  each  parameterization,  I  plot  the  Southern 
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Figure  3-2:  Steady  state  surface  Southern  Ocean  P04  sensitivity  to  global  dust  in¬ 
crease  (absolute  factor)  and  S.  Ocean  overturning  (Sv)  for  A)  net  scavenging  case,  B) 
scavenging-desorption  case,  and  C)  complexation  case.  For  the  net  scavenging  and 
scavenging-desorption  case,  an  increase  in  global  dust  supply  results  in  the  drawdown 
of  P04  with  little  sensitivity  to  the  strength  of  vertical  exchange.  For  the  complex¬ 
ation  case  (C),  P04  drawdown  is  muted,  with  greatest  drawdown  when  the  vertical 
exchange  flux  is  weakened  to  10  Sv. 


Ocean  surface  [P04]  as  a  function  of  a  global  increase  in  aeolian  iron  supply,  relative 
to  today’s,  and  for  several  rates  of  Southern  Ocean  vertical  mass  exchange.  By 
increasing  the  dust  flux  10  times  globally,  surface  [PO4]  is  depleted  in  both  the  net 
scavenging  and  scavenging-desorption  models  (Figure  3-2  A  and  B).  There  is  little 
sensitivity  to  the  strength  of  Southern  Ocean  overturning.  In  strong  contrast,  for  the 
complexation  parameterization  (Figure  3-2  C),  even  with  global  dust  increase  of  10 
times  and  the  strength  of  vertical  exchange  decreased  by  50%,  it  is  not  possible  to 
completely  drawTdown  surface  [PO4]  in  this  model. 

Comparing  the  importance  of  Fe  supplied  to  the  euphotic  zone  by  dust  to  upwelled 
Fe  gives  insight  into  the  underlying  mechanistic  differences.  I  plot  the  fraction  of  iron 
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directly  supplied  by  dust  to  the  surface  Southern  Ocean  (Figure  3-3)  and  the  deep 
water  dissolved  iron  concentration  (Figure  3-4)  for  the  three  models.  For  the  net 
scavenging  and  scavenging-desorption  parameterization,  the  fraction  of  iron  supplied 
to  the  Southern  Ocean  euphotic  zone  regionally  by  dust  is  small  (5-10%).  I  find 
that  this  fraction  does  not  respond  to  increased  global  aeolian  dust  supply  (Figure 
3-3  A  and  B)  because  the  slow  net  scavenging  rate  enables  iron  derived  from  low- 
latitude  dust  to  be  transported  at  depth  to  the  deep  Southern  Ocean.  Therefore 
the  upwelled  source  of  iron  from  the  Southern  Ocean  increases  in  proportion  to  the 
global  dust  deposition  (Figure  3-4  A  and  B),  tracking  the  total  dust  supply.  For 
the  complexation  parameterization,  the  fraction  of  iron  supplied  by  dust  increases 
strongly  with  aeolian  dust  deposition  (Figure  3-3  C).  This  is  because  the  imposed, 
finite  ligand  concentration  places  an  upper  limit  on  the  deep  water  iron  concentration 
(Figure  3-4  C)  and  therefore  on  the  upwelled  iron  source. 

It  is  possible  that  ligand  production  increases  as  a  function  of  increased  dust 
flux,  as  evidenced  by  Rue  and  Bruland  (1997)  during  the  Iron-Ex  II  study  in  the 
equatorial  Pacific.  As  my  sensitivity  study  using  the  complexation  model  with  an 
elevated  [Lj]~  4  nM  shows  (see  Chapter  2),  deep  water  [Fer]  would  increase  with 
total  ligand  concentration  and  so  might  the  upw’elling  supply.  However  I  have  not 
parameterized  this  specific  mechanism  here. 

The  three  different  parameterizations  of  deep  water  iron  cycling  are  able  to  capture 
the  observed  distribution  of  iron  in  the  modern  ocean.  However  these  parameteriza¬ 
tions  lead  to  very  different  sensitivities  of  surface  phosphate  drawdown  in  conditions 
of  increased  dust  supply.  It  is  perhaps  premature  to  suggest  that  one  parameterization 
is  more  realistic  than  another  in  this  regard,  though  the  complexation  parameteri¬ 
zation  resolves  more  details  of  the  system  as  it  is  presently  understood.  This  result 
is  very  significant  for  model  projections  of  glacial-interglacial  biogeochemical  change, 
such  as  that  of  Watson  et  al.  (2000)  which  applied  a  scavenging  based  parameter¬ 
ization.  Clearly  it  is  imperative  to  continue  to  seek  more  observational  data  and  a 
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Figure  3-3:  Percent  Iron  in  surface  Southern  Ocean  derived  from  dust 
(dust/dust+upwelling)  for  the  A)  net  scavenging  case,  B)  scavenging  and  desorption 
case  and  C)  complexation  case  as  a  function  of  dust  flux  and  S.  Ocean  overturning. 


Global  Ocean  Dust  Increase 


Figure  3-4:  Southern  Ocean  deep  FeT  response  to  global  ocean  dust  increase  and 
upwelling  strength  for  A)  net  scavenging  case,  B)  scavenging-desorption  case,  and 
C)  complexation  case.  For  the  net  scavenging  and  scavenging-desorption  case,  an 
increase  in  global  dust  supply  results  in  a  proportional  increase  in  deep  water  [Fex]- 
For  the  complexation  case  (C)  Southern  Ocean  deep  water  [Fex]  increase  is  much  less 
than  for  A  and  B. 
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deeper  understanding  of  the  key  processes  in  order  to  make  more  appropriate  models 
for  climate  change  studies. 
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Chapter  4 


Global  Iron  Cycling:  Simulations 
Using  an  Ocean  General 
Circulation  Model 

4.1  Introduction 

Explorations  with  a  multi-box  model  showed  that  three  different  parameterizations 
are  able  to  capture  the  observed  deep  water  Fe  gradients.  A  box  model  is  qualita¬ 
tively  informative  and  provides  a  simple  framework  to  focus  on  the  biogeochemical 
controls  on  iron  cycling.  Due  to  its  economical  efficiency,  a  box  model  also  allows  for 
the  exploration  of  parameter  values,  not  possible  with  a  three-dimensional  general 
circulation  model  because  of  computational  restraints.  But  there  are  limitations  to 
the  box  model,  as  it  does  not  resolve  intra-basinal  and  vertical  gradients.  In  addition, 
the  parameterization  of  physical  terms  such  as  advection  and  mixing  are  very  crude. 
In  order  to  gain  insight  into  the  role  of  physical  processes  on  iron  cycling,  as  well  as 
higher  resolution  results  than  box  model  simulations  can  provide,  I  implement  my  iron 
cycling  parameterizations  within  the  framework  of  a  three-dimensional  ocean  general 
circulation  model  (GCM).  Results  from  the  sensitivity  studies  with  the  box  model 
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guides  the  choice  of  parameter  values  in  the  general  circulation  model.  In  this  chapter, 
I  discuss  the  physical  properties  of  the  GCM,  the  biogeochemical  components  of  the 
model  and  the  aeolian  forcing  fields.  I  also  present  results  from  the  implementation  of 
the  A)  Net  Scavenging,  B)  Scavenging-Desorption  and  C)  Scavenging-Complexation 
iron  parameterizations  in  the  three-dimensional,  coarse  resolution  general  circulation 
model. 


4.2  Physical  Model 

The  MIT  GCM  is  configured  at  coarse  resolution  (2.8x2.8  degrees,  15  levels)  globally 
(Marshall  et  al.  1997a,  b;  Adcroft  et  al. ,  1997).  The  model  is  forced  with  clima¬ 
tological  surface  wind  stresses.  Surface  heat  flux  is  imposed  from  climatology  with 
an  additional  relaxation  toward  observed  SST,  and  surface  salinity  is  relaxed  toward 
climatology.  In  these  coarse  resolution,  global  ocean  circulation  studies,  mesoscale 
eddy  transfer  effects  are  achieved  using  schemes  related  to  the  parameterization  of 
Gent  and  McWilliams  (1990).  Vertical  turbulent  mixing  in  the  surface  mixed  layer 
of  the  ocean  is  parameterized  using  convective  adjustment  (Dutay  et  al,  2002). 

In  order  to  assess  the  physical  controls  on  Fe  cycling  in  the  global  ocean,  it  is 
necessary  that  the  the  general  circulation  model  adequately  represents  major  phys¬ 
ical  processes  occurring  in  the  ocean.  In  this  section,  I  briefly  describe  the  residual 
circulation  of  the  model  and  present  meridional  sections  of  temperature  and  salinity 
in  the  Atlantic  and  Pacific  basins. 

4.2.1  Residual  Circulation 

Modern  ocean  models  transport  tracers  using  the  ‘residual  circulation’  -  the  net  ef¬ 
fect  of  advection  by  mean  flow  and  a  contribution  from  induced  circulation  (Gent  and 
McWilliams,  1990;  Karsten  et  al.,  2002;  Marshall  et  al.,  2002).  The  residual  circu¬ 
lation  of  this  configuration  is  presented  in  Figure  4-1.  The  overturning  is  strongest 
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Figure  4-1:  The  Global  Residual  Mean  of  the  model.  Fluxes  are  in  Sv  (106  m3  s_1) 
Figure  courtesy  of  Takamitsu  Ito. 
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Figure  4-2:  Modeled  vertical  velocity  (m  yr_1)  at  depth  of  50  meters.  Upwelling  is 
occurring  in  yellow-shaded  areas.  Figure  courtesy  of  Mick  Follows. 


in  the  North  Atlantic  at  ~45°N  due  to  the  formation  of  North  Atlantic  Deep  Wa¬ 
ter.  The  residual  circulation  of  the  model  is  weaker  in  the  Southern  Ocean  than 
estimates  of  the  meridional  overturning  circulation  from  data  and  inverse  modeling 
suggest  (Karsten  and  Marshall,  2002;  Sloyan  et  al,  2001).  The  wind  driven  gyres  are 
apparent  in  the  upper  layers. 

Vertical  Velocity 

Upwelling  is  an  important  transport  pathway  for  nutrients  from  the  deep  to  return 
to  the  euphotic  zone.  Near  the  surface,  upwelling  is  strongly  related  to  wind  stress 
forcing  and  Ekman  pumping.  In  Figure  4-2, 1  plot  the  vertical  velocity  of  the  model  at 
50  meters.  As  seen  in  the  figure,  the  ocean  gyres  are  areas  of  downwelling,  suppressing 
the  return  of  nutrients  from  the  deep.  This  results  in  low  biological  productivity  in 
the  ocean  gyres.  Major  areas  of  upwelling  are  the  Southern  Ocean,  the  equatorial 
Pacific  and  the  subpolar  gyres. 

4.2.2  Temperature  and  Salinity  Sections 

Figure  4-3  illustrates  sections  of  modeled  temperature  in  the  Atlantic  (330°E)  and 
the  Pacific  (198°E)  compared  to  observed  temperature  sections  (Levitus  and  Boyer, 
1994).  In  both  basins,  the  model  captures  the  structure  of  the  thermocline.  In  the 
equatorial  regions  and  subtropics  of  both  basins,  the  modeled  thermocline  is  too 
shallow,  resulting  in  lower  temperatures  than  observed.  In  the  North  Atlantic,  the 
plume  of  anomalously  warm  water  at  55°N  suggests  that  deep  mixing  is  too  vigorous 
in  this  region  of  the  model. 

In  Figure  4-4,  I  present  sections  of  modeled  salinity  in  the  Atlantic  and  the  Pa¬ 
cific.  The  high  salinity  signature  of  North  Atlantic  Deep  Water  and  the  low  salinity 
signature  of  the  Antarctic  Intermediate  Water  are  clearly  visible,  although  modeled 
values  are  slightly  higher  than  observations  (Levitus  et  al,  1994).  The  model  agrees 
well  with  observations  in  the  subtropical  waters  of  the  Atlantic.  Although  modeled 
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(A)  Observed  Temperature:  Atlantic  (330E) 


(B)  Observed  Temperature:  Pacific  (198E) 
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(C)  Modeled  Temperature:Atlantic  (330  E) 
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(D)  Modeled  Temperature:  Pacific  (198E) 
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Figure  4-3:  Comparison  of  observed  and  model  temperature.  Meridional  section  of 
observed  temperature  (°C,  Levitus  and  Boyle,  1994)  in  the  A)  Atlantic  basin  at  330°E 
and  B)  Pacific  basin  at  199°E.  Model  temperature  for  the  C)  Atlantic  basin  and  D) 

the  Pacific  basin. 


salinity  is  higher  than  observed  throughout  the  Pacific,  the  model  is  able  to  capture 
the  structure  of  the  salinity  field. 


4.3  Biogeochemical  Component 


In  addition  to  iron,  the  biogeochemical  tracers  explicitly  carried  in  my  model  are 
phosphate  (P04  )  and  dissolved  organic  phosphorus  (DOP).  Biological  uptake  and 
regeneration  are  indexed  to  phosphorus.  The  governing  equations  for  P04  and  DOP 
in  the  biogeochemical  model  at  the  surface  are: 
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(A)  Observed  Salinity:  Atlantic  (330  E)  (B)  Observed  Salinity:  Pacific  (198  E) 
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(C)  Modeled  Salinity:Atlantic  (330°E)  (D)  Modeled  Salinity:  Pacific  (198°E) 
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Figure  4-4:  Comparison  of  observed  and  model  salinity.  Meridional  section  of  ob¬ 
served  salinity  (ppt,  Levitus  et  al.,  1994)  in  the  A)  Atlantic  basin  at  330°E  and  B) 
Pacific  basin  at  199°E.  Model  salinity  for  the  C)  Atlantic  basin  and  D)  the  Pacific 
basin. 
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dP04 

dt 

dDOP 

dt 


-  V  -(u P04)  +  v  •  («  V  PO4 )  -  r  +  A DOP 

-  V  -(u DOP)  +  V  *  (K  V  POP)  +  i/r  -  XDOP 


r  =  //P04 


/ 

+  F$ 1  +  I0 


(4.1) 

(4.2) 

(4.3) 


At  depth,  the  equations  are: 


dP^4  =  -  V  -(uP04)  +  v  •  («  V  PO4)  +  (1  -  ^)r  +  XDOP  -  — (4-4) 
=  v  •  (u DOP)  +  V  •  («  V  -DOP)  -  XDOP  (4.5) 

t/  6 

r° 

Export  —  J  (1  —  P)rdz  (4-6) 

J  hmix 

F(z)  =  Export(-  —-)~b  (4-7) 

hmix 

In  the  above  equations,  u  is  the  transformed  Eulerian  mean  velocity  and  k  is  a 
mixing  tensor  representing  isopycnal  mixing  following  Gent  and  McWilliams  (1990). 
A  third-order  upwind,  flux  corrected  advection  scheme  is  used  for  tracers  (Dutkiewicz 
et  ai,  2001). 

r  represents  the  biological  uptake,  which  is  limited  by  light,  phosphate,  and  iron. 
In  conditions  where  Fe  and  light  are  replete,  I  assume  surface  PO4  to  be  the  limiting 
nutrient  which  is  utilized  with  a  characteristic  timescale,  1/ p,  of  about  1  month.  Iron 
limitation  is  represented  by  Michaelis-Menten  kinetics.  The  half  saturation  constant 
for  iron  ( Ks )  is  globally  uniform  but  is  adjusted,  within  the  range  of  measured  values 
(Price  et  al.,  1994;  Fitzwater  et  ai,  1996),  to  optimize  the  modeled  surface  [PO4]  and 
[Fex]  distributions. 

Two-thirds  of  exported  nutrient  (u)  enters  the  surface  dissolved  organic  phos¬ 
phorus  (DOP)  pool.  The  imposed  timescale  for  remineralization  of  DOP  (1/A)  is  6 
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months.  One-third  of  biological  uptake  [(1-  u)  T]  is  rapidly  exported  as  particulate 
at  depth  to  the  PO4  pool  (Yamanaka  and  Tajika,  1997)  using  an  empirical  power 
law  relationship  determined  by  Martin  et  al.  (1987).  The  depth  of  the  euphotic  zone 
( hmix )  is  75  meters  and  the  power  law  coefficient  (6)  is  0.9. 

The  light  field  varies  as  a  function  of  time  and  latitude  by  the  following  equation: 


I  =  Q0{  1  -  o)  cos(z) 
cos(z) 


D,  .  cos  S cos  6 sin  D. 

— [sin  0  sin  (p-\ - 

7T  D 


D  =  cos  1  (—  tan  5  tan  <f>) 


(4.8) 

(4.9) 
(4.10) 


The  solar  constant  (Q0)  is  1,367  Wm“2.  The  albedo  (a)  varies  latitudinally  and 
temporally  according  to  North  et  al.  (1981).  D/tt  is  the  the  sunlit  region  as  a 
fraction  of  latitude  circle  and  the  term  in  brackets  in  equation  4.9  is  the  average 
cosine  zenith  angle  over  the  sunlit  region.  Solar  declination  ( S )  is  calculated  daily 
using  the  astronomical  formula  of  Paltridge  and  Platt  (1976).  The  half-saturation 
constant  for  light  ( I0 )  is  30  Wm'2. 


4.3.1  Aeolian  Flux  Forcing  Field 

The  primary  source  of  iron  to  the  global  ocean  is  aeolian.  Gao  et  al  (2001)  create 
seasonal  aeolian  iron  flux  maps  extrapolated  from  in  situ  marine  boundary  layer 
measurements.  Based  on  the  average  concentration  of  Fe  in  continental  crust,  they 
assume  3.5  weight  percent  of  dust  is  iron.  Dust  deposition  is  strongest  in  the  Northern 
Hemisphere  with  maximum  deposition  in  the  subtropical  and  tropical  Atlantic  and 
the  North  Pacific  (Figure  4-5).  The  Fe  flux  is  strongest  during  the  summer  months 
in  the  Atlantic.  The  trade  winds  carry  Saharan  dust  westward  depositing  Fe  across 
the  Atlantic.  In  the  North  Pacific,  dust  deposition  is  highest  during  the  spring.  The 
prevailing  westerlies  transport  dust  originating  from  the  Asian  desert  into  the  North 
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Pacific.  In  contrast,  the  eastern  Southern  Pacific  and  the  Southern  Ocean  receive 
very  little  dust  throughout  the  year. 

The  model  is  forced  seasonally  to  capture  the  spatial  and  temporal  variability  in 
dust  deposition.  Based  on  explorations  with  the  box  model  (see  Chapter  2),  I  assume 
1%  of  the  iron  entering  the  surface  ocean  is  soluble  (Jickells  and  Spokes, 2001). 


4.3.2  Iron  Parameterization 


The  three  iron  paramcterizations  described  and  tested  in  Chapter  2  using  a  multi-box 
model  framework  are  implemented  here  in  the  ocean  general  circulation  model  with 
the  exception  of  the  description  of  scavenging.  In  the  box  model  simulations,  scaveng¬ 
ing  was  modeled  as  a  first  order  process,  limited  only  by  [Fe/].  Laboratory  experiments 
using  thorium  and  oceanic  field  observations  indicate  that  particle  concentration  also 
limits  scavenging  (Honeyman  et  al.,  1998  and  references  therein).  Compiling  oceanic 
field  data,  Honeyman  et  al.  (1988)  fit  a  power  law  function  to  describe  the  relationship 
between  scavenging  rate  (ksc)  and  particle  concentration  (Cp): 

k$c  =  koCp  (4-11) 

where  k0  represents  the  scavenging  rate  when  particles  are  not  limiting,  Cp  is 
the  particle  concentration  and  <j>  is  the  slope.  I  employ  this  relationship  to  model 
scavenging  as  a  function  of  Cp.  In  the  model,  [Cp]  is  calculated  for  every  grid  point  for 
each  level  from  the  modeled  biological  flux.  Thus,  ksc  is  faster  in  biologically  active 
regions  of  the  ocean  and  since  [Cp]  decreases  with  depth,  the  modeled  scavenging 
rate  becomes  slower  at  depth.  Since  the  empirical  relationship  was  calculated  using 
thorium,  I  scale  the  scavenging  rate  for  iron  such  that  the  average  surface  scavenging 
rate  falls  within  the  range  used  in  the  box  model  simulations. 
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4.4  Results 


While  all  three  parameterizations  successfully  reproduced  the  broad  observed  pat¬ 
terns  of  deep  water  iron  gradients  in  the  box  model  simulations,  the  ocean  general 
circulation  model,  with  a  more  sophisticated  representation  of  physics,  identifies  some 
differences.  In  this  section,  I  present  results  for  the  three  different  iron  parameteri¬ 
zations  within  the  context  of  an  ocean  general  circulation  model:  A)  net  scavenging 
case,  B)  scavenging-desorption  case,  and  C)  scavenging-complexation  case. 

For  each  case,  I  compare  the  modeled  PO4  distributions  to  maps  of  observed  [PO4] 
(Figure  4-6,  Conkright  et  al,  1994).  Due  to  the  paucity  of  iron  measurements  in  the 
ocean,  it  is  not  possible  to  make  detailed  comparisons  between  modeled  [Fer]  and 
observations.  To  assess  the  model’s  ability  to  reproduce  observed  [Fe^]  distributions, 
I  compare  the  [Fe^]  with  the  compilation  of  Fe^  measurements  at  various  depths 
presented  (Figure  1-4)  in  Chapter  1. 

4.4.1  Net  Scavenging  Results 

In  the  context  of  the  box  model,  the  net  scavenging  case  was  the  simplest  description 
of  iron  that  resolved  the  deep  water  iron  gradients.  I  set  the  variable  scavenging 
scaling  parameter  (t )  to  0.0035  and  test  whether  this  simple  description  is  able  to  still 
capture  the  observed  global  distribution  within  the  context  of  a  more  sophisticated 
description  of  ocean  circulation  and  transport  processes. 

Phosphate 

Figure  4-7  shows  the  [P04]  distribution  at  the  surface,  935m,  and  2495m  using  the 
net  scavenging  parameterization.  In  the  Southern  Ocean  and  major  portions  of  the 
Pacific,  the  model  does  predict  excess  (i.e.  nutrient  concentration  greater  than  the 
limiting  nutrient  multiplied  by  the  Refield  ratio)  surface  [P04],  but  does  a  poor  job  of 
distinguishing  the  HNLC  equatorial  region  from  the  low  PO4  regions  of  the  subtrop- 
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ical  gyres.  This  may  be  due  to  a  lack  of  resolution  of  important  physical  processes 
(Aumont  et  al.,  1999).  Coarse  resolution  ocean  GCMs,  such  as  this,  do  not  resolve 
important  equatorial  dynamical  processes,  resulting  in  an  overestimation  of  upwelling 
velocity.  This  has  been  pointed  out  as  the  underlying  cause  of  ‘nutrient  trapping’  in 
such  models  with  surface  nutrient  restoring  (Najjar  et  al.,  1992).  Since  I  do  not  re¬ 
store  P04  to  observations  at  the  surface,  the  excess  P04  brought  to  the  surface  by 
upwelling  may  be  advected  laterally,  resulting  in  high  [PO4]  in  the  subtropics.  In 
addition,  due  to  the  high  fraction  of  biologically  utilized  PO4  transfered  to  DOP  and 
the  slow  remineralization  rate  of  DOP,  the  lifetime  of  DOP  is  very  long  in  the  upper 
layers  of  the  model.  This  can  result  in  the  remineralization  of  DOP  in  the  upwelling 
water,  near  to  the  surface,  resulting  in  high  [PO4]  in  surface  waters.  Additionally,  the 
flux  of  iron  to  the  subtropical  Pacific  may  be  higher  than  the  prescribed  dust  forc¬ 
ing  in  the  model.  This  is  further  discussed  in  Chapter  5.3.1.  While  the  magnitude 
of  excess  [PO4]  is  lower  than  observed  in  the  Northern  Pacific,  the  model  does  not 
completely  deplete  [PO4]  in  this  area.  At  depth,  the  model  reproduces  the  observed 
deep  water  [PO4]  gradients. 

Iron 

The  [Fey]  is  depleted  in  the  surface  ocean,  except  in  the  high  aeolian  flux  regions  - 
the  North  Atlantic,  the  North  Pacific  and  Indian  Ocean.  The  concentration  in  the 
North  Atlantic  is  much  higher  than  observations  (Figure  4-8).  The  model  may  be 
missing  processes  that  occur  in  the  surface  ocean.  This  is  discussed  in  more  detail  in 
section  4.4.3.  While  the  model  successfully  reproduces  the  observed  deep  water  [Fey] 
gradients  at  935  meters  and  2495  meters,  the  [FeT]  is  much  higher  than  observed  in 
the  Atlantic.  In  the  Indo-Pacific  basin  and  the  Southern  Ocean,  the  [Fer]  is  in  general 
agreement  with  the  observations,  but  higher  by  ~0.2  nM. 

The  high  [Fe^]  signature  of  the  Atlantic  water  is  carried  by  North  Atlantic  Deep 
Water  into  the  the  Indian  sector  of  the  Southern  Ocean.  There  are  no  observations 
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Figure  4-7:  Modeled  [P04]  (ju) M  at  A)  surface,  B)  935m  and  C)  2495m  for  the  net 
scavenging  case. 
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from  this  area  at  depths  greater  than  1000m  against  which  to  compare  the  model 
results.  This  high  [Fe^]  at  depth  is  much  higher  than  observations.  Allowing  the  net 
scavenging  rate  to  increase  would  improve  the  deep  water  [Fer]  results  in  the  North 
Atlantic,  but  it  would  also  result  in  lowering  the  [Fer]  in  the  other  basins,  which  are 
consistent  with  observations. 

While  the  net  scavenging  case  does  reproduce  the  observed  deep  water  [Fer]  gra¬ 
dients,  it  is  not  able  to  deplete  Fe  in  the  Atlantic  waters  within  the  context  of  this 
more  sophisticated  physical  framework.  Perhaps  complexation  and/or  desorption  are 
necessary  to  accurately  describe  the  global  iron  cycle. 


4.4.2  Scavenging-Desorption  Results 

In  the  box  model,  the  scavenging-desorption  parameterization  was  able  to  reproduce 
the  broad  features  of  the  large  scale  distribution  of  dissolved  iron.  Due  to  technical 
difficulties  with  the  sinking  flux,  I  was  not  able  to  test  whether  this  parameterization 
is  still  able  to  capture  the  deep  water  iron  gradients  within  a  context  of  a  more 
sophisticated  circulation.  This  is  work  I  plan  to  pursue  in  the  future. 


4.4.3  Complexation  Results 

Of  the  three  parameterizations,  the  complexation  model  is  the  most  detailed  and 
explicitly  describes  Fe-ligand  interactions.  Within  the  context  of  the  box  model  sim¬ 
ulations,  it  is  the  only  parameterization  that  explicitly  accounts  for  iron’s  speciation. 
In  the  GCM  simulation,  I  set  the  variable  scavenging  scaling  parameter  (r)  to  0.125 
for  the  complexation  case  and  ligand  strength  to  log(A')  =  11  based  on  the  results  of 
sensitivity  studies  using  the  box  model  presented  in  chapter  2. 
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Phosphate 


Figure  4-9  illustrates  the  modeled  phosphate  distribution  compared  to  global  PO4 
maps  (Conkright  et  al,  1994),  Figure  4-6).  The  model  captures  the  broad  pattern 
of  the  observed  phosphate  distribution,  most  notably  excess  [P04]  in  the  Southern 
Ocean.  Although  the  excess  [PO4]  is  lower  than  observed  in  the  Southern  Ocean, 
this  is  an  improvement  over  earlier  nutrient  cycling  models  that  did  not  include  iron. 
It  was  necessary  to  retain  surface  PO4  to  observed  values  (Orr,  2002)  or  set  a  long 
export  timescale  in  order  to  return  excess  PO4  in  the  Southern  Ocean  (Me  Kinley  et 
al,  2000).  Here  it  is  explicitly  achieved.  In  the  South  Pacific,  the  model  is  not  able 
to  capture  the  PO4  surface  gradients  between  the  equatorial  and  subtropical  regions. 
Possible  explanations  are  discussed  in  section  4.4.1.  In  the  thermocline  and  below 
2000m,  the  model  successfully  reproduces  the  magnitude  and  direction  of  the  deep 
water  gradients. 

Iron 

In  qualitative  agreement  with  observations,  the  model  predicts  elevated  surface  [For] 
in  the  Atlantic  basin  (Figure  4-10).  This  is  the  case  because  the  flux  of  iron  to  the 
ocean  is  highest  in  the  Atlantic.  The  magnitude  of  the  modeled  [Fer]  in  the  Atlantic 
is  much  higher  than  observed  and  exceeds  the  solubility  of  iron  (Liu  and  Millero, 
2002),  which  should  result  in  precipitation  of  iron.  Also,  formation  of  colloids  and 
aggregation  of  iron  may  be  occurring  in  surface  waters  (Wu  et  al,  2001).  These 
processes  are  not  included  here  and  I  suggest,  may  explain  the  offset  in  the  magnitude 
of  excess  [Fer]  between  observations  and  modeled  results.  I  plan  to  include  them  in 
future  work. 

Additionally,  some  studies  have  suggested  (Fridlind  and  Jacobsen,  2000;  Siefert, 
1999;  Willey  et  al,  2000)  that  the  further  iron  travels  in  the  atmosphere,  the  more  it 
is  chemically  processed,  increasing  its  solubility  when  deposited  in  the  ocean.  Thus 
Fe  deposited  near  its  source  may  be  less  soluble  than  iron  that  has  traveled  farther. 
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Figure  4-9:  Modeled  [PO4]  (/xM)  at  A)  surface,  B)  935m  and  C)  2495m  for  the 
complexation  case. 
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Figure  4-10:  Modeled  [Fer]  (nM)  at  A)  surface,  B)  935m  and  C)  2495m  for  the 
complexation  case. 
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The  solubility  of  the  dust  may  also  be  affected  by  the  iron  concentration  in  surface 
waters.  In  iron  rich  waters,  the  ability  to  dissolve  the  aeolian  derived  iron  may  be 
lower  if  the  water  is  already  close  to  saturation. 

The  model  also  predicts  excess  surface  iron  in  the  North  Pacific  and  Indian  Ocean, 
reflecting  areas  of  elevated  dust  flux  (see  Figure  4-5),  in  good  agreement  with  obser¬ 
vations.  In  the  rest  of  the  Indo-Pacific  basin  and  the  Southern  Ocean,  [Fex]  is  close 
to  complete  drawdown  as  is  to  be  expected  since  these  are  the  traditional  high  nu¬ 
trient,  low  chlorophyll  regions  of  the  ocean.  Broadly,  the  modeled  surface  patterns 
are  generally  consistent  with  observations,  but  quantitative  differences  do  exist,  most 
notably  high  [Fex]  in  the  Atlantic.  Not  allowing  [Fex]  to  exceed  solubility  limits  and 
considering  the  role  of  colloids  and  aggregation  would  likely  improve  model  results  in 
surface  waters  and  are  mechanistically  reasonable.  I  plan  to  study  controls  on  surface 
[Fex]  in  the  future. 

Figure  4-10B  compares  the  modeled  [Fer]  distribution  at  935  meters  with  mea¬ 
surements.  In  good  agreement  with  the  observations,  [Fex]  is  highest  in  the  North 
Atlantic  and  North  Pacific.  Modeled  [Fex]  is  also  high  in  the  Northern  Indian  Ocean, 
but  there  are,  as  yet,  no  observations  to  compare  with.  The  model  successfully  re¬ 
produces  the  pattern  of  decreasing  [Fex]  in  the  Pacific  moving  from  north  to  south 
and  predicts  lowest  [Fex]  in  the  Southern  Ocean,  in  agreement  with  the  few  available 
measurements.  There  are  no  published  measurements  of  iron  in  the  South  Pacific 
and  South  Atlantic,  making  it  impossible  to  assess  the  performance  of  the  model  for 
these  regions. 

There  are  very  few  Fe  measurements  at  depths  greater  than  2000  meters  (see  Fig¬ 
ure  1-4,  bottom  panel).  Observations  suggest  highest  concentrations  in  the  Atlantic 
basin  and  the  North  Pacific.  The  rest  of  the  Pacific  and  the  Southern  Ocean  appear 
to  have  lower  values,  ranging  between  0.4-0.6  nM.  My  model  results  show  good  agree¬ 
ment  with  the  few  measurements  that  exist.  The  model  predicts  a  maximum  [Fer] 
in  the  North  Atlantic  of  0.7  nM,  reflecting  the  high  aeolian  input.  The  signature 
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of  iron  carried  by  North  Atlantic  Deep  Water  is  visible  in  the  Atlantic  sector  of  the 
Southern  Ocean.  Since  dust  deposition  is  very  low  in  the  Southern  Ocean,  advection 
appears  to  play  an  important  role  in  delivering  iron  to  the  deeper  Southern  Ocean 
waters.  In  the  model,  deep  water  [Fey]  appear  to  be  uniform,  whereas  observations 
suggest  a  decreasing  gradient  from  north  to  south.  More  observational  data  is  needed 
to  constrain  the  model  at  depth,  especially  in  the  South  Pacific  and  Southern  Ocean. 

Unlike  the  other  two  parameterizations,  the  complexation  case  is  also  able  to 
predict  the  speciation  of  the  Fe.  Figure  4-11  shows  the  ratio  of  FeL  to  Fey  in  the 
model.  In  good  agreement  with  observations,  over  90%  of  the  Fe  is  organically  bound. 
Only  in  the  North  Atlantic,  where  surface  [Fey]  are  higher  than  the  imposed  [Ly]  (1 
nM),  is  the  FeL/Fey  ratio  low. 

Comparison  to  previously  published  models 

I  compare  my  results  to  Archer  and  Johnson’s  (2000)  case  with  a  two-ligand  iron 
model  to  a  3-dimensional  ocean  general  circulation  model  (described  in  Chapter  2). 
Their  modeled  surface  PO4  field  (Figure  4-12)  for  the  two  ligand  case  does  a  better  job 
of  predicting  elevated  [PO4]  in  the  Southern  Ocean.  However,  in  order  to  achieve  this 
result,  Archer  and  Johnson  (2000)  have  decreased  the  solubility  of  the  dust  reaching 
the  Southern  Ocean  relative  to  the  rest  of  the  ocean.  Yet,  preliminary  evidence 
(Fridlind  and  Jacobsen,  2000;  Siefert,  1999;  Willey  et  al.,  2000)  suggest  the  longer 
dust  remains  in  the  atmosphere,  the  more  soluble  it  becomes  due  to  the  low  pH  of 
clouds.  This  then  implies  that  dust  reaching  the  Southern  Ocean  would  have  a  higher 
solubility. 

At  2500  meters  depth,  the  [Fey]  is  uniform  in  Archer  and  Johnson’s  (2000)  two 
ligand  case  (Figure  4-13).  The  deep  water  [Fey]  is  ~0.55  nM.  In  the  southeastern 
coast  off  the  coast  of  South  America  the  Archer  and  Johnson  (2000)  model  predicts 
elevated  [Fey].  My  model  better  captures  the  gradient  between  the  Atlantic  and 
Pacific  basin.  At  the  time  of  Archer  and  Johnson’s  (2000)  work,  it  was  thought  that 
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Figure  4-12:  Modeled  surface  PO4  distribution  for  Archer  and  Johnson’s  (2000)  two- 
ligand  case  assuming  6%  of  Mahowald  et  al.  (1999)  is  soluble  (Figure  taken  from 
Archer  and  Johnson,  2000). 

deep  water  [Fey]  were  the  same  in  both  basins.  Thus,  Archer  and  Johnson  (2000) 
tuned  their  model  with  this  information  in  mind.  Newer  measurements  show  that  a 
gradient  does  exist  between  the  two  basins. 

While  my  model  has  a  similar  iron  parameterization,  differences  in  the  solution 
can  be  attributed  to  the  values  chosen  for  certain  parameters.  Archer  and  Johnson’s 
(2000)  strong  ligand  has  a  very  high  stability  constant  (log(/\  )=13)  and  a  concentra¬ 
tion  of  0.6  nM  ensuring  a  saturated  ligand.  Additionally,  Archer  and  Johnson  (2000) 
use  a  scavenging  rate  of  0.8  yr-1,  four  times  greater  than  the  scavenging  rate  I  use. 
With  such  a  fast  scavenging  rate,  any  iron  that  is  not  complexed  is  quickly  removed 
from  the  system.  Since  [Lx]  =  0.6  nM,  this  essentially  sets  the  iron  concentration  to 
be  uniform,  but  is  no  longer  consistent  with  the  data.  I  suggest  that  Archer  and  John¬ 
son’s  (2000)  results  are  essentially  a  limit  case.  I  find  a  better  solution  if  the  strong 
constraints  of  high  log(A^)  and  low  [Lx]  are  relaxed.  Archer  and  Johnson  (2000)  also 
remove  ‘excess’  iron  (i.e.  iron  above  the  amount  needed  to  support  biological  pro¬ 
duction)  from  the  surface,  not  allowing  any  of  it  to  be  transported  at  depth.  While 
my  results  suggest  there  are  mechanisms  occurring  in  the  surface  waters  to  remove 
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Figure  4-13:  Modeled  Fex  distribution  at  350m  and  2500m  from  Archer  and  Johnson’s 
(2000)  two-ligand  model  case.  (Figure  taken  from  Archer  and  Johnson,  2000). 

excess  iron  such  as  precipitation  or  aggregation  of  colloidal  iron,  it  is  still  possible 
that  a  portion  of  the  excess  iron  is  transported  at  depth  with  the  sinking  of  the  North 
Atlantic  Deep  Water. 

4.5  Summary 

Using  a  coarse-resolution  ocean  general  circulation  model,  I  find  it  is  not  simple  to  re¬ 
produce  observed  Fe  gradients  and  structure  with  the  net  scavenging  and  scavenging- 
desorption  descriptions.  In  contrast,  the  complexation  model  easily  captures  the 
observed  Fe  distribution  without  much  tuning.  I  compare  the  results  of  the  com¬ 
plexation  case  with  those  of  Archer  and  Johnson  (2000).  I  suggest  that  their  model 
is  a  special  limit  case  and  that  by  relaxing  constraints,  my  model  is  better  able  to 
predict  deep  water  Fer  gradients.  Differences  in  the  solution  of  the  two  models  can  be 
attributed  to  how  each  model  accounts  for  excess  Fex  at  the  surface,  ligand  strength, 
and  scavenging  rate. 
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Chapter  5 


Regional  and  Global  Iron 
Distributions 


In  the  previous  chapter,  I  showed  that  the  complexation  parameterization  of  iron  is 
best  able  to  explain  global  iron  distributions.  In  this  chapter  I  further  analyze  the 
results  from  the  complexation  parameterization  of  iron.  I  define  a  tracer,  Fe*  that 
calculates  to  what  degree  a  water  mass  is  iron  limited.  I  also  calculate  the  residence 
time  of  iron  in  the  global  ocean  and  evaluate  the  importance  of  upwelling  as  a  source 
of  iron  to  the  euphotic  zone. 

In  contrast  to  the  last  chapter,  where  I  compared  my  model  results  to  observations 
globally,  I  assess  how  well  my  model  results  capture  temporal  and  regional  variations 
in  iron  distribution.  Lastly,  I  discuss  the  results  of  two  sensitivity  studies:  the  response 
of  modeled  global  FeT  distribution  to  A)  a  different  aeolian  forcing  field  and  B) 
parameterizing  the  uptake  of  Fe  and  P  as  a  ratio  that  varies  in  relation  to  ambient 


5.1  A  Tracer  of  Iron  Limitation:  Fe* 


The  scavenging  of  iron  from  the  water  column  results  in  the  decoupling  of  Fe  and 
PO4 .  Regions  where  upwelled  waters  are  highly  deficient  in  Fe  require  a  high  aeolian 
Fe  flux  to  compensate,  otherwise  iron  limitation  occurs.  Here,  I  construct  a  tracer 
which  tracks  the  relative  magnitude  of  decoupling  between  Fe  and  P04.  Similar 
to  N*  (Gruber  and  Sarmiento,  1997),  I  subtract  the  contribution  of  the  soft  tissue 
pump  from  the  Fe^  distribution  to  reveal  the  balance  between  physical  transport  and 
scavenging  of  iron.  For  the  surface,  I  define: 


Fe*  —  Fer  —  FtFeP04 


dFe * 
dt 

dFe * 
dt 


_  dFeT  n  dP04 
dt  Rpe  dt 

=  ocFin  -  u  ■  v(Fe*)  -  kscFe 


(5.1) 

(5.2) 

(5.3) 


At  depth: 

dFe* 

-ft-  =  -«  •  V{Fe*)  -  kscFe  (5.4) 

By  removing  the  soft  tissue  pump  contribution,  Fe*  reveals  the  balance  between 
advected  and  scavenged  Fe.  It  is  important  to  note  that  RFe,  the  ratio  of  Fe:P04 
must  be  specified  to  calculate  Fe*.  In  nature,  the  Fe:C  is  not  likely  to  be  constant 
(Sunda  and  Huntsman,  1995),  but  in  the  model,  I  have  assumed  a  constant  Fe:C  ratio 
of  4  pmokl  mol  and  a  Redfield  ratio  of  C:P  of  117:1,  resulting  in  an  Fe:P  ratio  of 
.47  mmolrl  mol.  A  positive  Fe*  implies  there  is  adequate  Fe  to  support  the  complete 
biological  utilization  of  P04,  while  a  negative  Fe*  means  that  there  is  a  deficit  in  Fe. 

In  Figures  5-1,  5-2  and  5-3, 1  have  plotted  zonally  averaged  sections  of  Fe*  for  the 
Atlantic,  Indo-Pacific  and  Southern  Ocean  basins,  respectively.  In  the  North  Atlantic, 
where  aeolian  input  is  high,  Fe*  is  above  zero  for  all  depths  (Figure  5-1).  The  signature 
of  the  strong  aeolian  deposition  is  carried  to  the  deep  waters  by  the  western  boundary 
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Figure  5-1:  Zonally  averaged  section  of  Fe*  in  the  Atlantic  basin. 

current  and  southward  by  North  Atlantic  Deep  Water  to  approximately  35°S.  Only 
southward  of  35°S  does  scavenging  deplete  the  excess  Fe  that  originated  from  N. 
Atlantic  surface  waters. 

Only  in  the  surface  waters  of  the  Indo-Pacific  basin  from  approximately  25°-40°N 
is  Fe*  above  zero  (Figure  5-2).  The  waters  upwelled  from  below  are  deficient  in  Fe 
(negative  Fe*),  but  aeolian  deposition  is  able  to  compensate  so  macronutrients  are  still 
strongly  drawn  down  in  the  surface.  Fe*  is  most  negative  in  the  deep  Northern  Pacific 
waters.  Since  these  are  the  oldest  waters  in  the  ocean,  the  iron  has  been  stripped  away 
by  scavenging,  while  PO4  has  been  accumulating,  resulting  in  the  largest  decoupling 
between  P04  and  Fe  in  the  ocean.  Therefore,  although  dust  deposition  is  high  in  the 
North  Pacific,  HNLC  conditions  are  still  prevalent.  In  the  Southern  Ocean  (Figure 
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Figure  5-2:  Zonally  averaged  section  of  Fc*  in  the  Indo-Pacific  basin 


Fe*:  Southern  Ocean 


Figure  5-3:  Zonally  averaged  section  of  Fe*  in  the  Southern  Ocean  basin 


5-3),  Fe*  is  negative  throughout  the  basin  and  decreases  with  depth,  reflecting  the 
dominance  of  scavenging  over  transport. 


5.2  Iron  Residence  Time 

For  geochemical  problems,  we  sometimes  consider  the  ocean  to  be  a  well  mixed  reser¬ 
voir  in  which  the  resulting  concentration  of  an  element  is  dependent  on  the  rate  of 
input  of  the  element  into  the  ocean  and  the  rate  of  reactions  that  remove  the  element 
from  the  ocean.  Assuming  steady  state,  the  input  (source)  and  output  (sink)  should 
balance.  Dividing  the  rate  of  input  or  output  by  the  amount  of  the  element  in  the 
ocean  provides  us  with  the  residence  time  (t),  or  mean  lifetime  of  that  element  in 
the  ocean,  which  is  also  an  indicator  of  the  amount  of  time  it  takes  for  an  element  to 
respond  to  a  change  in  its  input  or  output  rate. 

I  calculate  the  residence  time  of  Fe  in  the  model  and  discuss  how  it  compares  to 
other  estimates.  The  soluble  dust  flux  entering  the  ocean  is  ~  2*109  mol  Fe  yr_1  and 
the  amount  of  Fe  in  the  ocean  at  steady  state  averaged  for  one  year  is  5.76*  10n  mol 
Fe,  resulting  in  a  r  =  285  years.  Such  a  long  residence  time  implies  that  Fe  can  be 
transported  quite  far  prior  to  scavenging,  hence  identified  with  water  mass  tracers. 

Bacon  and  Anderson  (1982)  show  that  the  residence  time  of  a  metal  can  be  calcu¬ 
lated  by  considering  the  amount  of  iron  that  is  lost  due  to  scavenging.  They  consider 
the  ratio  of  the  adsorbed  metal  concentration  ([Me/>])  to  the  total  metal  concentra¬ 
tion  ([Mer])  to  be  an  estimator  of  residence  time.  Using  this  method,  Bruland  et 
al.  (1994)  estimate  a  residence  time  ranging  between  70-140  years  for  Fe.  While  my 
estimate  predicts  a  longer  residence  time  than  Bruland  et  al.  (1994)  for  Fe,  due  to 
the  large  uncertainties,  they  are  probably  not  inconsistent.  Bruland  et  al.  (1994) 
estimate  the  [Mcp]  by  measuring  the  fraction  of  Fe  that  is  leachable  by  acetic  acid. 
Since  the  pH  of  acetic  acid  is  lower  than  pH  of  the  sea,  it  may  leach  off  more  Fe  than 
was  actually  scavenged.  This  would  imply  that  Bruland  et  a  Vs  (1994)  estimate  is  a 
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lower  bound  estimate.  Additionally,  Bruland  et  al.  (1994)  make  a  global  residence 
time  estimate  using  samples  from  the  central  North  Pacific  gyre.  In  my  calculation, 
there  are  also  uncertainties  associated  with  the  aeolian  flux  of  iron  into  ocean.  Due 
to  the  episodic  nature  of  aeolian  Fe  input,  it  is  difficult  to  accurately  estimate  the 
global  average  yearly  dust  flux.  Furthermore,  the  solubility  of  aeolian  iron  is  very 
uncertain  in  seawater.  If  less  than  1%  of  aeolian  Fe  is  soluble,  as  preliminary  dust 
leaching  experiments  suggest  (Jickells  and  Spokes,  2001),  and  the  global  dust  flux 
I  used  is  an  overestimate,  my  estimate  can  be  viewed  as  an  upper  bound.  Hence, 
within  uncertainties,  the  estimates  are  consistent. 

5.3  Importance  of  upwelled  iron  to  the  euphotic 
zone 

The  two  sources  of  ’new’  iron  to  the  euphotic  zone  are  aeolian  deposition  and  up- 
welling.  The  amount  of  soluble  Fe  entering  the  ocean  in  the  model  is  ~2*109  mol 
Fe  yr-1.  The  amount  of  Fe  upwelled  into  the  euphotic  zone  is  1.3*109  mol  Fe  yr~\ 
accounting  for  40%  of  the  supply  to  the  surface.  My  calculation  agrees  well  with  the 
estimate  of  Fung  et  al.  (2000)  based  on  the  results  from  a  mixed  layer  model.  Archer 
and  Johnson  (2000),  though,  predict  that  ~70-75  %  of  iron  is  supplied  by  upwelling. 
Archer  and  Johnson’s  model  predicts  a  uniform  [FeT]  of  ~0.6nM  throughout  the  deep 
waters,  while  more  recent  observations  show  that  [FeT]  are  the  lower  in  the  Pacific 
and  Southern  Ocean.  Thus,  Archer  and  Johnson  (2000)  may  have  overestimated  the 
contribution  of  upwelling  as  a  source  of  iron  to  the  euphotic  zone. 
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Figure  5-4:  Modeled  profile  at  BATS  (35°N,  58°W)  (diamonds)  compared  to  observa¬ 
tions  (*)  for  the  month  of  July.  Observational  data  provided  by  Boyle  and  co-workers 
(unpublished) . 

5.4  Temporal  and  Vertical  Variations:  Compar¬ 
isons  to  Data 

In  this  section,  I  compare  my  model  results  with  observational  data  from  specific  sta¬ 
tions.  Since  the  deposition  of  Fe  varies  episodically  and  the  model  is  forced  seasonally, 
I  also  examine  the  fluctuation  of  [Fe^]  seasonally  in  the  surface  layer. 

Figure  5-4  displays  the  vertical  profile  of  ’dissolved’  (<  0.4/zm)  [Fe^]  from  35°N, 
58° W  (Wu  et  al,  2001),  as  well  as  the  modeled  profile.  The  modeled  [Per]  are  much 
higher  than  observed,  suggesting  either  that  the  soluble  aeolian  Fe  estimates  for  the 
Atlantic  are  too  high  or  there  are  unrepresented  sinks  in  the  model  for  iron  in  surface 
waters,  such  as  precipitation  or  aggregation  and  subsequent  sinking  of  Fe.  I  believe 
the  latter  is  the  more  likely  possibility.  Below  800m,  the  model  results  agree  well 
with  observations.  The  iron  profile  is  straight  below  1000m,  although  the  observed 
[Fey]  decreases  slightly  at  ~3700m.  The  maximum  [Fe^]  for  both  the  model  and 
observations  is  ~0.7nM. 
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Figure  5-5:  Modeled  profile  at  Station  ALOHA  (22°N,  158°W)  compared  to  observa- 
tions.  Observational  data  provided  by  Boyle  and  co-workers  (unpublished). 

At  station  ALOHA,  the  model  predicts  a  HNLC  region  since  [Fex]  is  completely 
depleted,  but  [PO4]  is  elevated.  In  contrast,  the  observations  show  excess  [Fex]  and 
depleted  [P04]  in  surface  waters  here  (Figure  5-5).  This  suggests  that  the  aeolian  forc¬ 
ing  used  in  the  model  is  too  weak  at  station  ALOHA  , ‘nutrient  trapping’  is  occurring, 
or  the  DOF  remineralization  timescale  is  too  long.  The  latter  two  possibilities  are 
discussed  in  further  detail  in  Section  4.4.1.  Below  300m,  the  model  shows  better  agree¬ 
ment  with  the  observations.  The  observations  show  a  maximum  [Fex]  at  ~1200m, 
which  varies  seasonally.  The  model  maximum  [Fex]  is  at  a  slightly  shallower  depth 
and  does  not  predict  [Fex]  as  high  as  observations  indicate.  This  discrepancy  can 
be  attributed  to  the  presence  of  the  Loihi  seamount,  a  local  source  of  hydrothermal 
iron.  At  3000m,  the  model  predicts  an  [Fex]  of  0.5  nM,  while  Boyle  (unpublished) 
measured  0.45  nM. 

The  modeled  profile  in  the  North  Pacific  (50°N,  145°W)  matches  the  observations 
from  Johnson  et  al.  (1997)  well  (Figure  5-6).  The  Fex  is  drawn  down  at  the  surface 
and  has  a  maximum  [Fex]  at  ~1000  meters.  At  intermediate  depths,  the  modeled  Fex 
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Figure  5-6:  Modeled  profile  in  the  North  Pacific  at  50°N,  145°W  (diamonds)  for  July 
compared  to  observations.  Observational  data  (*)  from  Johnson  et  al.  (1997). 

is  more  strongly  scavenged  than  the  observations.  Iron  from  the  continental  margins 
may  add  iron  at  depth  in  this  region,  possibly  explaining  why  the  modeled  [Fey]  is 
lower  at  intermediate  depths.  In  the  equatorial  Pacific,  the  model  completely  draws 
down  [Fey],  also  seen  in  the  observations  (Figure  5-7).  The  model  has  a  maximum 
[Fey]  at  1000m  due  to  remineralization,  and  decreases  at  intermediate  depths  due  to 
scavenging.  The  [Fey]  increases  at  depths  greater  than  3000m  because  the  scavenging 
rate  decreases  due  to  the  decrease  in  particle  density. 

In  Figure  5-8,  I  plot  the  modeled  profile  [Fey]  (diamonds)  at  56°S,  143 °E  in  the 
Southern  Ocean  and  compare  it  to  measurements  made  by  Sohrin  et  al.  (2000).  The 
model  depletes  Fe  completely  at  the  surface.  Measurements  by  Sohrin  et  al.  (2000) 
show  excess  Fe  at  the  surface,  but  Coale  et  al.  (2002)  have  shown  Fe  to  be  depleted  in 
surface  Southern  Ocean  waters.  At  intermediate  depths  (1000-2500m),  the  modeled 
concentration  is  slightly  higher,  but  has  a  similar  profile  to  the  measurements.  My 
model  successfully  predicts  a  lower  [Fey]  at  1000m  than  is  observed  in  the  Pacific  and 
Atlantic  basins.  Below  3000m,  the  model  agrees  well  with  observations. 
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Figure  5-7:  Modeled  profile  in  the  Equatorial  Pacific  at  3°S,  140°W  (diamonds) 
compared  to  observations  (*)  in  November.  Observational  data  from  Johnson  et 
al.  (1997). 
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Figure  5-8:  Modeled  profile  in  the  Southern  Ocean  at  56°S,  143°E  (diamonds)  com¬ 
pared  to  observations  (*).  Observational  data  from  Sohrin  et  al.  (2000). 
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In  summary,  beneath  the  surface  layer,  the  model  agrees  well  with  observed  profiles 
for  different  regions  of  the  ocean.  At  the  surface,  the  model  overestimates  the  excess 
[Fer]  in  high  flux  regions.  In  areas  with  lower  aeolian  input  such  as  the  Southern 
Ocean  and  the  central  Pacific,  the  model  is  able  to  drawdown  [Fe^]  completely,  while 
observations  do  show  a  low  [FeT]  at  the  surface.  In  the  equatorial  and  North  Pacific, 
both  HNLC  regions,  the  model  achieves  better  agreement  at  the  surface  with  the 
data.  Processes  missing  from  this  model  which  might  account  for  these  discrepancies 
include  a  variable  Fe:C  ratio,  or  a  more  complex  parameterization  of  biological  export 
that  accounts  for  the  varying  needs  of  different  class-size  phytoplankton.  Inclusion  of 
these  processes  may  improve  the  model’s  performance  at  the  surface  and  are  areas  I 
hope  to  explore  in  the  future. 

Surface  Transects 

In  the  western  North  Atlantic,  Boyle  et  al.  (unpublished)  collected  surface  samples 
along  45°W  between  KPN  and  35°N  in  January  2001  and  then  again  in  July  2001  to 
study  interannual  variability  in  surface  [Fey].  They  found  (Figure  5-9)  little  variation 
in  the  [Fey] -  In  January,  the  observations  show  [Fe^]  =  0.6  nM,  while  it  is  0.55 
nM  in  July  at  10°N.  The  model  predicts  a  [Fer]  =  0.78  nM  in  January  (Figure  5-9) 
and  0.82  nM  in  July  at  10°N.  Although  the  model  predicted  [Fer]  are  higher  than 
observations,  the  fact  that  [Fer]  at  the  surface  in  January  and  July  are  comparable, 
is  consistent  with  observations.  Modeled  [FeT]  are  much  higher  than  observations 
due  to  the  high  aeolian  flux  forcing  for  this  region.  Of  interest  is  the  surface  [FeT] 
between  January  and  July  (Figure  5-10).  Model  results  suggest  a  strong  seasonal 
cycle.  Monthly  measurements  of  surface  Fer  would  confirm  or  refute  this  suggestion. 

The  observations  and  modeled  results  are  somewhat  surprising,  since  aeolian  flux 
is  low  during  January,  but  high  in  July  (Figure  5-11).  The  dust  flux  estimates  used 
in  the  model  simulation  are  high  for  June  through  November,  but  drop  significantly. 
I  calculate  a  local  residence  time  of  35  days  for  iron  in  surface  waters  at  this  location. 
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Figure  5-10:  Monthly  modeled  surface  [Fe^]  at  10°N,  45°W. 

This  implies  that  the  soluble  Fe  in  the  dust  deposited  in  November  still  influences 
surface  waters  in  January. 

For  station  ALOHA,  Boyle  (unpublished)  compiled  seasonal  surface  [FeT]  and 
found  a  sharp  increase  in  surface  [Fey]  in  April  (Figure  5-12A).  Although  quanti¬ 
tatively  the  model  predicts  much  lower  surface  [FeT]  than  observed  (Figure  5-12B), 
the  model  is  able  to  reproduce  the  seasonal  cycle,  which  is  controlled  by  the  aeolian 
dust  flux,  with  highest  flux  during  the  spring  (Figure  5-13).  The  model  also  predicts 
high  [PO4]  for  station  ALOHA,  which  is  at  odds  with  observations  (Figure  4-6).  This 
suggests  the  aeolian  flux  used  to  force  the  model  is  too  low,  creating  a  HNLC  region 
where  there  is  not  one.  Alternatively  it  could  be  related  to  ’nutrient  trapping’  or  the 
long  timescale  for  remineralization  of  DOP  in  the  model.  Both  of  these  possibilities 
are  discussed  in  greater  detail  in  Section  4.4.1. 

The  model  does  a  poor  job  of  predicting  the  magnitude  of  [Fer]  for  the  surface 
transects  I  compared  my  model  to.  In  the  Atlantic,  the  model  values  are  higher  than 
observed  and  in  the  Pacific,  the  model  predicts  lower  [Fer]  than  observed.  Improved 
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Figure  5-11:  Monthly  dust  flux  at  10°N,  45°W  used  in  the  model. 

dust  flux  estimates  and  including  surface  processes  such  as  precipitation  and  aggrega¬ 
tion  of  Fe  may  improve  model  results.  However,  the  aim  of  the  model  was  to  capture 
the  cycling  of  Fe  in  the  deep  water,  so  it  is  not  entirely  surprising  that  the  model  can 
improve  in  the  surface.  The  model  is  able  to  successfully  mimic  the  seasonal  cycle  at 
Station  ALOHA  and  suggests  there  is  a  strong  seasonal  surface  signal  in  the  western 
Atlantic  also. 

5.5  Sensitivity  to  Aeolian  Forcing 

Aeolian  deposition  is  the  primary  source  of  iron  into  the  global  ocean  and  thus  plays 
an  important  role  in  setting  the  ocean  Fe  distribution.  Yet  estimates  of  aeolian  Fe 
deposition  are  poorly  constrained  due  to  the  episodic  nature  of  dust  deposition.  I  test 
the  sensitivity  of  the  ocean  iron  distributions  to  the  range  of  estimates  in  regional 
dust  flux  by  forcing  the  model  with  a  second  set  of  aeolian  flux  forcing  fields  that 
was  generated  from  an  atmospheric  dust  transport  model  (Mahowald  et  al,  2003). 
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Figure  5-13:  Monthly  dust  flux  at  Station  ALOHA  used  in  the  model. 

In  this  section,  I  describe  the  aeolian  flux  field  and  discuss  the  corresponding  results 
of  my  model. 

Forcing  Field 

Mahowald  et  al.  (2003)  generated  daily  global  aeolian  flux  maps  from  a  22-year 
atmospheric  model  simulation.  Figure  5-14  displays  the  seasonal  global  iron  flux 
used  in  this  sensitivity  study.  I  assume  3.5  weight%  of  dust  is  iron  and  1%  solubility. 
Similar  to  the  Gao  forcing  field,  aeolian  iron  deposition  is  highest  in  the  Northern 
Atlantic.  Deposition  is  also  high  in  the  Arabian  Sea  and  the  Northwestern  Pacific, 
although  relative  to  Gao  et  al.  (2001)  (see  Figure  4-5),  the  Northwestern  Pacific 
flux  is  weaker  in  the  Mahowald  et  al.  (2003)  simulation.  The  most  notable  difference 
between  the  two  aeolian  flux  forcing  fields  is  the  strength  of  deposition  in  the  Southern 
Hemisphere.  The  Mahowald  et  al.  (2003)  simulation  suggests  much  higher  dust 
deposition  in  the  Southern  Hemisphere,  especially  off  of  the  Australian  coast  and  the 
northeastern  coast  of  South  America. 
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Figure  5-14:  Seasonal  dust  deposition  (mg  Fe  m  2  yr  *)  simulated  from  an  atmo¬ 
spheric  model  (Mahowald  et  al.,  2003). 
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Results 


The  Southern  Ocean  surface  [P04]  is  lower  (Figure  5-15)  when  the  model  is  forced 
with  the  Mahowald  et  al.  (2003)  fields,  attributable  to  the  higher  aeolian  flux  in  the 
Southern  Hemisphere,  relative  to  the  Gao  et  al.  (2001)  forcing. 

At  the  surface,  the  Mahowald  simulation  is  able  to  drawdown  FeT  throughout  the 
Pacific  basin  and  the  Southern  Ocean  (Figure  5-16).  The  magnitude  and  extent  of 
excess  Fer  are  higher  for  the  Mahowald  simulation  in  the  North  Atlantic  and  the 
Indian  Ocean,  while  the  band  of  elevated  surface  [Fer]  in  the  North  Pacific  produced 
by  the  Gao  simulation  is  not  present  in  the  Mahowald  results.  Therefore,  the  [P04] 
is  also  higher  in  the  surface  North  Pacific  of  the  Mahowald  simulation. 

At  1000m,  both  simulations  predict  a  decreasing  [Fer]  gradient  from  the  North 
Pacific  to  the  central  Pacific,  in  broad  agreement  with  the  observations.  The  North 
Pacific  [Fer]  is  higher  in  the  Gao  simulation  because  there  is  a  stronger  aeolian  input 
coming  from  Asia  than  in  the  Mahowald  simulation. 

At  2500  meters,  the  solution  to  both  simulations  looks  remarkably  similar.  El¬ 
evated  Fer  transported  by  the  North  Atlantic  Deep  Water  is  carried  further  in  the 
Mahowald  simulation,  likely  due  to  the  higher  aeolian  input  into  this  region  in  the 
Mahowald  forcing  field.  Significantly,  the  variations  in  the  two  aeolian  flux  fields  are 
not  large  enough  to  affect  the  model  solution  in  the  deep  water  [Fer]  •  While  qualita¬ 
tively,  the  two  model  results  are  similar  for  the  surface  [P04]  and  [Fer]  distributions, 
the  magnitudes  differ  considerably.  The  Gao  simulation  predicts  higher  surface  [PO4] 
in  the  Southern  Ocean,  in  better  agreement  with  observations  (although  still  lower 
than  observations)  than  the  Mahowald  simulation.  The  Mahowald  simulation  results 
for  surface  [PO4]  in  the  equatorial  and  subtropical  Pacific  are  closer  to  observations. 
The  P04  distribution  resulting  from  Mahowald  simulation  in  the  Southern  Hemi¬ 
sphere  reflects  the  observed  gradient  between  the  subtropical  and  subpolar  waters. 
The  simulation  is  not  able  to  differentiate  these  water  masses  by  their  [P04].  Both 
simulations  overestimate  excess  surface  iron  concentrations  in  the  North  Atlantic  and 
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Indian  Ocean,  very  likely  reflecting  the  absence  of  processes  that  control  [Fe?-]  in 
surface  waters. 

It  is  important  not  to  interpret  either  forcing  field  as  better  than  the  other.  I  could 
have  tuned  the  solubility  parameter  (a)  to  improve  the  model  output,  but  chose  not 
to  in  order  to  understand  the  sensitivity  of  [Fer]  to  temporal  and  spatial  differences 
in  aeolian  input.  For  example,  decreasing  ct,  the  solubility  parameter,  would  improve 
the  Southern  Ocean  surface  [PO4]  in  the  Mahowald  case,  but  it  would  also  increase 
the  surface  P04  throughout  the  surface  waters  of  the  model  domain,  making  the 
model  results  look  particularly  unrealistic  in  the  equatorial  and  subtropical  Pacific 
waters. 


5.6  Sensitivity  to  Variable  Fe:C  Uptake  Ratio 

In  order  to  optimize  growth,  Hudson  and  Morel  (1990)  suggest  phytoplankton  growing 
in  low  Fc  environments  must  increase  their  surface  arearvolume  ratio  or  decrease  their 
cellular  iron  requirement  to  optimize  growth.  Oceanic  environments  that  are  severely 
iron  limited,  such  as  the  Southern  Ocean,  are  dominated  by  small  cells  (Miller  et 
ai,  1991),  supporting  Hudson  and  Morel’s  (1990)  conclusion  with  respect  to  cell 
size.  Sunda  and  Huntsman  (1995)  found  lower  cellular  iron  requirements  for  oceanic 
phytoplankton  compared  to  coastal  phytoplankton,  supporting  Hudson  and  Morel’s 
(1990)  second  conclusion.  In  Figure  5-17, 1  plot  the  results  of  Sunda  and  Huntsman’s 
(1995)  experiment  for  oceanic  phytoplankton  at  ambient  [Fer]  <  15  nM.  The  log-log 
plot  shows  a  weak  correlation  of  increasing  Fe:C  with  increasing  [Fer].  I  fit  a  power 
law  to  the  data: 


R  =  R0[Fe]n  (5.5) 

where  log(i?o)  (the  Y-intercept)  equals  1.18  and  n,  the  slope  equals  0.44.  It  is  impor¬ 
tant  to  note  that  the  lowest  ambient  [Fer]  in  Sunda  and  Huntsman’s  (1995)  exper- 
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Figure  5-16:  Modeled  [FeT]  (nM)  using  the  dust  forcing  from  Mahowald  et  al.  (2003) 
at  (A)  surface,  (B)  935m  and  (C)  2495m. 


Fe:C  vs.  [Fey]  (Sunda/Huntsman,  1995  (intracellular,  Fe<15nM) 


Figure  5-17:  Variance  in  intracellular  Fe:C  ratio  (log(Fe:C)*10-6)  as  a  function  of 
ambient  [Fey]  (nM).  The  relationship  is  very  weak  and  the  experiment  was  conducted 
at  high  ambient  [Fex]  (Sunda  and  Huntsman,  1995). 

iments  were  2  nM.  In  the  open  ocean,  rarely  have  [Fey]  above  1  nM  been  observed. 
Nevertheless,  I  impose  this  relationship  in  my  model,  allowing  Fe:C  ratio  to  vary 
based  on  the  ambient  [Fex]-  I  assume  phytoplankton  reach  a  point  where  it  is  not 
possible  to  increase  the  Fe:C  efficiency  (arbitrarily  set  to  1*10“6  mol  Fe:mol  C  ).  The 
P:C  ratio  is  constant  (117:1).  In  the  reference  case,  the  Fe:C  ratio  is  4*10-6  mol  Fe:l 
mol  C. 


Results 

In  Figure  5-18,  I  plot  the  [P04]  at  the  surface,  935  m,  and  2495  m.  By  allowing 
the  Fe:C  to  decrease  with  decreasing  surface  [Fer]  in  accordance  with  Sunda  and 
Huntsman’s  (1995)  data,  the  Pacific  switches  from  Fe  limitation  to  P04  limitation. 
While  the  model  is  still  Fe  limited  in  the  Southern  Ocean  and  the  equatorial  Pacific, 
the  excess  [PO4]  is  much  lower  than  observed.  Due  to  higher  export  production, 


110 


modeled  [PO4]  is  higher  than  observed  at  depth  in  the  old  Pacific  waters. 

Using  a  variable  Fe:C  ratio  has  only  a  small  effect  on  the  Fe  distribution  (Figure 
5-19).  [Fex]  is  lower  in  the  North  Atlantic,  but  elevated  for  the  rest  of  the  surface 
ocean.  In  my  reference  case,  I  use  an  Fe:C  ratio  of  4*10~6:1.  Using  the  power  law 
derived  from  Sunda  and  Huntsman’s  (1995)  work,  the  Fe:C  ratio  is  4*1CT6  only  when 
[Fer]  exceeds  2  nM.  Modeled  surface  [Fer]  only  reaches  such  concentrations  in  the 
North  Atlantic  surface  waters  and  the  Indian  Ocean.  Thus,  in  the  N.  Atlantic  and 
the  Indian  Ocean,  the  Fe:C  ratio  is  higher  in  the  variable  Fe:C  simulation  relative 
to  the  constant  Fe:C  case.  For  the  other  areas  of  the  ocean,  the  reverse  is  true.  At 
depth,  the  structure  is  essentially  unchanged,  but  the  [Fer]  is  slightly  lower  at  depth. 
This  may  be  attributed  to  less  FeT  being  biologically  utilized  at  the  surface  and  then 
subsequently  exported  at  depth. 

Simply  using  a  variable  Fe:C  ratio  in  the  model  did  not  improve  the  model  results. 
However,  other  free  parameters  could  be  tuned  in  concert.  Simultaneously  decreasing 
the  solubility  of  the  aeolian  iron  entering  the  ocean  may  have  improved  the  model’s 
response.  Leaching  studies  of  Fe  from  dust  do  suggest  that  less  than  1%  (the  solubility 
of  dust  in  my  model)  of  the  Fe  may  be  soluble  (Jickells  and  Spokes,  2001;  Spokes  and 
Jickells,  1996 ).  Different  size  class  species  may  have  different  Fe:C  ratios.  Rather  than 
parameterizing  export  production  simply,  an  ecosystem  model  that  accounts  for  the 
different  requirements  of  various  sized  phytoplankton  may  be  able  to  better  predict 
surface  [Fer]  and  [PO4].  While  experimental  evidence  does  suggest  Fe:C  ratios  are 
dependent  on  ambient  surface  water  [Fer]  (Sunda  and  Huntsman,  1995),  the  nature 
of  this  relationship  is  not  clear. 


5.7  Chapter  Summary 

I  examine  more  closely  the  complexation  case  by  comparing  my  model  results  to  ob¬ 
served  profiles  and  surface  transects  to  assess  the  ability  of  the  model  to  reproduce 
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Figure  5-18:  Modeled  [P04]  (/jM)  using  a  variable  Fe:P  ratio  based  on  the  data  from 
Sunda  and  Huntsman  (1995)  at  (A)  surface,  (B)  935m  and  (C)  2495m. 
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Figure  5-19:  Modeled  [Fey]  (nM)  using  a  variable  Fe:P  ratio  based  on  the  data  from 
Sunda  and  Huntsman  (1995)  at  (A)  surface,  (B)  935m  and  (C)  2495m. 


vertical  gradients  and  seasonal  variability.  Below  300m,  the  model  can  successfully 
predict  [Fer],  At  the  surface,  the  model  predicts  seasonal  variability  in  [Fer].  It  over¬ 
estimates  surface  [Fer]  in  high  aeolian  flux  regions  such  as  the  North  Atlantic.  This 
could  be  attributed  to  the  absence  of  processes  such  as  aggregation  and  subsequent 
sinking  or  precipitation  that  may  serve  as  sinks  for  iron.  The  model  also  predicts  iron 
limitation  for  the  subtropical  Pacific,  not  consistent  with  observations.  The  model 
may  be  underestimating  the  dust  flux  to  this  region,  misrepresenting  the  dynamics 
of  the  region  due  to  the  coarse  resolution  of  the  grid,  or  trapping  nutrients  beneath 
the  mixed  layer  as  discussed  in  Section  4.4.1. 

I  define  a  tracer,  Fe*  that  quantifies  the  degree  to  which  a  water  mass  is  iron 
limited.  At  all  depths,  the  North  Atlantic  has  a  positive  Fe*  value,  indicating  it  is 
never  iron  limited.  The  excess  iron  at  the  surface  due  to  aeolian  transport  is  carried 
southward  at  depth  by  North  Atlantic  Deep  Water.  Except  for  surface  waters  between 
25°N  and  40°N,  Fe*  is  negative  throughout  the  Indo-Pacific  basin.  The  positive  Fe* 
at  the  surface  in  the  North  Pacific  can  be  attributed  to  the  elevated  aeolian  flux 
coming  from  Asia.  The  deep  North  Pacific  has  the  most  negative  Fe*  values,  reflecting 
the  stripping  out  of  Fe  from  these  old  waters.  Southern  Ocean  waters  are  negative 
everywhere,  reflecting  the  low  aeolian  input  in  the  surface  waters,  and  the  dominance 
of  scavenging  over  transport  at  depth. 

Lastly,  I  calculate  the  residence  time  and  estimate  the  contribution  of  Fe^  to  the 
surface  waters  from  upwelling  and  aeolian  input.  I  estimate  the  residence  time  of  Fe 
to  be  ~285  years,  confirming  that  transport  plays  an  important  role  in  controlling 
deep  water  [Fer].  Globally,  upwelling  accounts  for  40%  of  ‘new’  iron  reaching  the 
euphotic  zone. 
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Chapter  6 


GCM  simulation:  Southern  Ocean 
Phosphate  Sensitivity  to  Increased 
Aeolian  Dust  Flux 


I  force  my  complexation  iron  parameterization  with  estimated  global  dust  fluxes  from 
the  Last  Glacial  Maximum  (LGM)  (Mahowald  et  al,  1999)  to  study  the  response  of 
Southern  Ocean  surface  [P04].  In  Chapter  3,  I  carried  out  the  same  simulation  using 
the  multi-box  model.  I  repeat  this  simulation  in  the  more  sophisticated  physical 
description  of  the  general  circulation  model  to  investigate  the  intrabasinal  response 
of  P04  in  the  Southern  Ocean.  Due  to  varying  wind  stresses  and  fluxes  of  salt  and 
heat  to  the  ocean,  as  well  as  extent  of  ice  cover  in  the  Southern  Ocean,  it  is  likely  that 
the  strength  of  upwelling  and  thermohaline  circulation  differed  considerably  during 
the  LGM.  I  do  not  consider  such  changes  here,  but  hope  to  explore  the  sensitivity 
of  Southern  Ocean  surface  PO4  drawdown  to  circulation  changes  in  the  future.  The 
box  model  results  provide  me  with  an  integrated  answer  and  are  not  able  to  resolve 
variations  spatially.  In  this  chapter  I  present  the  results  from  the  GCM  simulation. 
I  also  discuss  the  importance  of  upwelling  relative  to  aeolian  flux  as  a  source  of  new 
Fer  to  the  euphotic  zone  and  the  distribution  of  Fe*. 
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6.1  Aeolian  Forcing  Field 


Mahowald  et  al.  (1999)  used  linked  terrestrial  biosphere,  dust  source  and  atmospheric 
transport  models  to  simulate  the  dust  cycle  in  the  atmosphere  for  the  Last  Glacial 
Maximum.  Figure  6-1A  shows  the  Fe  forcing  field  used  in  the  increased  dust  flux 
scenario  (referred  to  as  ‘paleo’  simulation)  and  compare  it  to  the  forcing  field  used 
for  the  ‘modern’  simulations  (Figure  6-1B)  .  As  in  the  ‘modern’  simulation,  I  assume 
1%  of  aeolian  Fe  is  soluble  and  3.5  weight  percent  of  aeolian  dust  is  Fe.  The  model  is 
forced  seasonally.  In  comparison  to  the  modern  aeolian  iron  forcing  field,  the  ocean 
global  input  is  nine  times  higher  for  the  ‘paleo’  simulation  and  63  times  higher  in  the 
Southern  Ocean  (Fig  6-1C).  The  increase  in  dust  flux  is  within  the  range  of  increase 
that  ice  cores  and  sediments  suggest  (Petit  et  al,  1999;  Rea,  1994),  but  at  the  high 
end  of  estimates.  In  the  Northern  Hemisphere,  largest  increases  arc  (Figure  6-1C)  in 
the  North  Pacific  and  North  Atlantic.  In  the  Southern  Hemisphere,  the  dust  flux  is 
much  stronger  off  of  the  eastern  coasts  of  Australia  and  South  America. 

6.2  Global  Distributions 

In  this  section,  I  present  surface  P04  and  deepwater  FcT  distributions  forcing  the 
model  with  the  ‘paleo’  dust  fluxes  and  compare  them  to  the  solution  using  ‘modern’ 
dust  fluxes.  As  in  the  increased  dust  flux  simulation  using  the  box  model,  the  [Lr]  is 
constant  and  set  to  InM. 

6.2.1  Phosphate 

Surface  phosphate  drawdown  is  complete  over  large  regions  of  the  ocean,  yet  in 
the  high  nutrient,  low  chlorophyll  regions  of  the  ocean,  there  are  pockets  of  excess 
[P04](Fig  6-2A).  Intermediate  and  deep  water  [P04]  (Figure  6-2B,C)  are  higher  than 
the  model  predicted  for  the  ‘modern’  simulation,  due  to  the  export  of  P04  biologically 
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Figure  6-1:  A)  Aeolian  Flux  forcing  for  the  increased  dust  flux  simulation  (Mahowald 
et  al,  1999).  B)  The  dust  forcing  used  for  the  ‘modern’  simulations  presented  in 
Chapters  4  and  5  (log  scale).  C)  The  difference  between  between  the  ‘paleo’  and 
‘modern’  forcing  fields  (absolute  scale). 
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utilized  at  the  surface.  This  results  in  the  redistribution  of  surface  P04  at  depth. 


6.2.2  Iron 

Iron  is  completely  drawn  down  only  in  the  equatorial  Pacific  and  the  Southern  Ocean 
(Figure  6-3 A),  two  areas  that  are  Fe  limited  today.  Surface  [Fer]  is  above  5  nM  in 
certain  areas  of  the  Indian  Ocean  and  North  Atlantic.  Unless  [LT]  increases  as  a 
function  of  ambient  [Fex],  the  surface  [Fex]  exceeds  solubility  limits  (Liu  and  Millero, 
2002).  As  discussed  in  Chapter  4,  processes  that  may  be  occurring  which  lower  surface 
[Fex]  are  not  included  in  the  model. 

At  intermediate  depth  (Figure  6-3B),  the  deepwater  [Fex]  increases  in  all  basins. 
The  imposed  [Lx]  places  an  upper  limit  on  deep  water  [Fex].  Only  in  the  North 
Atlantic,  where  dust  input  is  highest  does  [Fer]  reach  1  nM.  In  the  North  Pacific  at 
935  meters,  [Fex]  increases  by  0.2  nM,  while  in  the  Southern  Ocean,  increases  as  high 
as  0.4  nM  are  observed.  North  Atlantic  Deep  Water  transports  Fe  into  the  Atlantic 
and  Indian  sectors  of  the  Southern  Ocean  (Figure  6-3C).  In  the  North  Pacific,  [Fer] 
only  increases  by  M).05  nM  relative  to  the  ‘modern’  simulation.  Since  the  deep  waters 
of  the  North  Pacific  are  the  oldest,  it  has  been  subjected  to  scavenging  longer  than 
any  other  water  mass. 


6.3  Southern  Ocean  Results 

Due  to  the  Southern  Ocean’s  potential  ability  to  absorb  pC02, 1  examine  more  closely 
the  effect  of  increased  dust  flux  on  surface  Southern  Ocean  [P04].  Specifically,  I 
analyze  the  change  in  Fe*  and  calculate  the  importance  of  upwelling  relative  to  aeolian 
flux  as  a  source  of  ‘new’  Fe  to  the  euphotic  zone. 
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6.3.1  Comparison  of  ‘modern’  and  ‘paleo’  surface  phosphate 
distribution 

The  average  Southern  Ocean  surface  [PO4]  decreases  from  —0.6  /xM  during  the  ‘mod¬ 
ern’  simulation  to  —0.1  /xM  for  the  ‘paleo’  simulation  (Figure  6-4A,  B).  The  box 
model  results  predicted  a  A  PO4  of  —0.25  /xM  (Figure  3-2C)  with  a  10  fold  increase 
in  global  dust  supply.  In  the  box  model  simulation,  dust  flux  was  increased  10  fold 
in  each  basin,  whereas  in  the  GCM-paleo  simulation,  the  increase  in  dust  flux  to  the 
Southern  Ocean  is  63  times  larger,  while  in  the  central  Pacific  dust  flux  only  increased 
—3  times.  Additionally  the  physics  in  the  box  model  are  very  idealized.  Due  to  these 
differences,  the  box  model  and  GCM  results  are  not  inconsistent. 

Although  PO4  is  completely  drawn  down  in  portions  of  the  Southern  Ocean  for 
the  ‘paleo’  simulation,  it  is  important  to  note  in  the  ‘modern’  simulation  that  model 
surface  [P04]  is  lower  than  observed.  The  observed  average  [P04]  in  the  surface 
Southern  Ocean  is  —2  /xM,  but  the  model  only  predicts  an  excess  [P04]  =  0.6  /xM. 
Thus,  the  results  should  not  be  interpreted  in  terms  of  the  absolute  [PO4]  for  the 
paleo  simulation,  rather  according  to  the  A  PO4.  Archer  and  Johnson  (2000)  use  the 
Mahowald  et  al.  (1999)  ‘paleo’  dust  forcing  and  find  that  surface  P04  only  decreases 
1012  moles  in  the  surface,  while  my  model  predicts  a  —2.5  times  greater  decrease  in 
total  surface  P04.  The  discrepancy  relates  to  the  large  uncertainty  in  dust  fluxes. 
Compared  to  the  Gao  et  al.  (2001)  ‘modern’  dust  flux,  the  Mahowald  ‘paleo’  fluxes 
are  8.75  times  higher  averaged  globally,  whereas  compared  to  the  Mahowald  et  al. 
(1999)  ‘modern’  fluxes,  the  Mahowald  ‘paleo’  fluxes  are  only  1.75  times  higher.  On 
the  other  hand,  Archer  and  Johnson  (2000)  carried  out  another  simulation  in  which 
atmospheric  dust  flux  was  increased  uniformly  over  the  ocean  surface.  A  four  fold 
increase  in  dust  flux  resulted  in  the  surface  inventory  of  PO4  decreasing  by  —3  /xM 
and  a  decrease  in  pC02  on  the  order  of  20-30  ppm,  closer  in  agreement  with  my 
results. 
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Figure  6-5:  Zonally  averaged  Fe*  in  the  Southern  Ocean  for  the  ‘paleo’  case. 

6.3.2  Fe* 

Fe*  changes  considerably  when  the  global  dust  flux  is  increased.  In  Figure  6-5,  I 
plot  the  zonally  averaged  calculated  Fe*  for  the  Southern  Ocean.  In  the  present  day, 
Fe*  is  negative  throughout  the  Southern  Ocean  (Figure  5-3).  Forcing  the  model  with 
LGM  aeolian  flux  forcing  fields  results  in  a  positive  Fe*  in  the  Southern  Ocean  surface 
waters,  attributed  to  the  increased  aeolian  input,  as  well  as  the  upwelling  of  high  Fe* 
North  Atlantic  Deep  Waters  to  form  Antarctic  Intermediate  Waters.  Although  Fe*  is 
positive  everywhere  in  the  surface,  PO4  is  not  completely  drawn  down.  This  implies 
that  light  is  also  limiting  productivity.  Negative  Fe*  values  are  found  only  in  the 
Antarctic  Bottom  Waters.  This  shift  from  negative  Fe*  waters  to  positive  Fe*  waters 
throughout  the  Southern  Ocean  suggests  that  the  balance  between  scavenging  and 
transport  has  shifted  in  favor  of  transport. 
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6.3.3  Local  Source  of  Iron  to  the  Southern  Ocean 

For  the  ‘modern’  results,  I  calculate  that  upwelling  supplied  4.65*  107  mol  Fe  yr-1  to 
the  euphotic  zone  in  the  Southern  Ocean.  Forcing  the  model  with  the  ‘paleo’  aeolian 
iron  fields  results  in  an  order  of  magnitude  increase  of  upwelled  iron  reaching  the 
euphotic  zone.  Due  to  the  increase  in  aeolian  flux  though,  upwelling  only  supplies 
22%  of  ‘new’  iron  reaching  the  euphotic  zone,  whereas  in  the  ‘modern’  simulation, 
upwelling  accounted  for  47  %  of  the  iron  reaching  the  euphotic  zone.  Evidence  sug¬ 
gests  (COHMAP  Members,  1988)  that  during  the  LGM,  the  wind  strength  increased, 
which  would  result  stronger  patterns  of  upwelling  and  downwelling  associated  with 
Ekman  processes,  and  hence  communication  between  surface  and  deep  waters.  On  the 
other  hand,  ice  cover  also  increased,  preventing  waters  from  upwelling  to  the  surface. 
My  results  do  not  take  either  of  these  effects  into  account.  Using  a  two-dimensional 
model  of  the  Southern  Ocean  developed  by  Taka  Ito,  I  plan  to  explore  the  effects  of 
such  changes  on  P04  drawdown.  Nevertheless,  my  results  suggest  that  increased  dust 
flux  changes  the  balance  between  upwelling  and  aeolian  flux  as  sources  of  ‘new’  iron 
to  the  euphotic  zone. 

6.4  Comment 

There  has  been  considerable  focus  on  seeding  the  Southern  Ocean  (Lam  and  Chisholm, 
2002)  with  iron  to  absorb  anthropogenic  CO2,  yet  preliminary  results  from  this  work 
and  others  (Archer  and  Johnson,  2000;  Lefevre  and  Watson,  1999;Watson  et  al. , 
2000)  suggest  that  increased  efficiency  of  the  biological  pump  due  to  iron  loading 
might  account  for  only  25-40  ppm  of  the  80  ppm  drawdown  of  glacial-interglacial 
change.  This  would  result  in  the  sequestration  of  ~100  gigatons  of  carbon,  equivalent 
to  only  10-15  years  of  current  anthropogenic  CO2  emissions. 
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Chapter  7 
Conclusions 


This  study  focused  on  understanding  the  controls  of  deep  water  iron  distribution 
in  order  to  assess  the  decoupling  of  iron  and  phosphorus  in  the  ocean  interior  and 
implications  for  upwelling  as  a  source  of  new  ‘iron’  to  the  euphotic  zone.  In  the 
beginning  of  the  thesis,  I  present  three  candidate  parameterizations  of  deep  ocean 
iron  cycling.  All  three  parameterizations  are  able  to  reproduce  the  observed  deep 
water  iron  gradients  in  the  context  of  an  idealized  six  box  model.  Using  a  3-D  ocean 
general  circulation  model,  which  has  a  more  sophisticated  representation  of  physics, 
identifies  some  differences.  The  model  supports  the  notion  that  the  binding  of  Fe  to 
an  organic  ligand  is  able  to  counteract  the  loss  of  Fe  due  to  scavenging,  allowing  [FeT] 
to  exceed  solubility  limits  (Liu  and  Millero,  2002)  and  broadly  captures  a  number  of 
significant  observable  details  including  deep  water  iron  gradients,  and  excess  ligand 
concentrations. 

I  define  a  tracer,  Fe*  that  indicates  the  degree  to  which  a  water  mass  is  iron 
limited  relative  to  P04.  Due  to  the  high  aeolian  flux  and  sinking  of  North  Atlantic 
Deep  Water,  the  Atlantic  basin  has  a  positive  Fe*  value  at  all  depths,  indicating  it 
is  never  Fe  limited.  HNLC  areas  are  also  regions  of  upwelling,  so  deep  water  Fe  is 
an  important  source  to  the  euphotic  zone.  In  the  HNLC  areas,  Fe*  is  negative  below 
the  euphotic  zone,  reflecting  the  dominance  of  scavenging  over  transport  at  depth 
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which  leads  to  the  decoupling  between  Fe  and  P04.  At  the  surface  in  HNLC  regions, 
the  aeolian  Fe  flux  is  not  large  enough  to  compensate  for  the  upw'elled  negative  Fe* 
waters. 

I  estimate  that  the  global  ocean  residence  time  of  iron  is  ~285  years.  Such  a  long 
residence  time  confirms  that  transport  plays  an  important  role  in  controlling  deep 
water  [Fer]  and  explains  why  modeled  deep  water  [Fer]  in  the  Southern  Ocean  is 
higher  than  can  be  supported  by  the  aeolian  flux  to  the  region. 


There  is  much  interest  in  assessing  the  importance  of  upwelling  as  a  source  of 
’new’  iron  to  the  euphotic  zone.  Having  developed  a  model  that  is  broadly  consistent 
with  the  limited  observations  of  deep  water  Fe  observations,  I  am  able  to  estimate 
the  contribution  of  iron  by  upwelling  to  the  euphotic  zone.  I  find  that  in  the  model 
40%  of  ’new’  iron  is  supplied  to  the  surface  waters  by  upwelling  in  good  agreement 
with  Fung  et  al.  (2000),  but  low  compared  to  Archer  and  Johnson  (2000).  This 
discrepancy  may  be  explained  by  the  higher  than  observed  deep  water  [Fey]  Archer 
and  Johnson’s  (2000)  model  predicts  for  the  Pacific  basin  and  the  Southern  Ocean. 

Due  to  the  potential  ability  of  iron  to  change  the  efficiency  of  the  carbon  pump  in 
the  remote  Southern  Ocean,  I  study  the  response  of  Southern  Ocean  surface  P04  to 
increased  aeolian  dust  flux.  My  box  model  results  suggest  that  a  global  ten  fold  in¬ 
crease  in  dust  flux  can  support  a  P04  drawdown  of  ~0.25/xM,  while  the  GCM  results 
suggests  a  P04  drawdown  of  0.5//M  in  the  Southern  Ocean.  Given  the  approxima¬ 
tion  and  uncertainties  in  each  model,  this  difference  is  perhaps  not  surprising.  The 
resolution  of  vertical  structure  in  the  GCM  may  be  a  significant  difference.  The  in¬ 
terpretation  of  these  results  with  regards  to  pC02  drawdown  should  be  made  with 
caution,  due  to  uncertainties  to  parameters  such  as  wind  strength,  extent  of  ice  cover 
and  the  response  of  ligand  production  to  increased  aeolian  iron  input. 
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7.1 


Next  Steps 


7.1.1  Research  Recommendations  from  a  Modeling  Perspec¬ 
tive 

While  the  complexation-scavenging  parameterization  of  iron  is  able  to  capture  the 
observed  deep  water  gradients,  comparisons  with  surface  iron  observations  show  that 
the  model  does  a  poor  job  of  predicting  surface  [Fer]-  These  results  suggest  that 
colloidal  chemistry  and  precipitation  may  be  added  to  the  model.  Research  efforts 
aimed  at  studying  the  formation  characterization,  and  cycling  of  Fe  colloids  would  be 
immensely  helpful. 

The  aeolian  flux  of  iron  is  poorly  constrained.  Global  monitoring  stations  and 
developing  the  capability  to  place  iron  sensors  on  buoys  in  remote  waters  would  aid 
in  better  quantifying  the  magnitude  and  periodicity  of  aeolian  iron  input.  In  addition, 
the  chemistry  of  aeolian  iron  in  clouds  and  its  solubility  in  marine  surface  waters  is 
poorly  constrained.  A  better  understanding  of  these  processes  would  be  useful  for 
such  models. 

Sunda  and  Huntsman  (1995)  suggest  that  marine  phytoplankton  adjust  their  Fe 
requirements  based  on  ambient  iron  concentrations.  Since  the  exact  nature  of  this 
relationship  is  not  known,  it  is  very  difficult  to  parameterize  this  effect,  as  seen  in 
the  sensitivity  study  in  Chapter  4.  Culture  work  aimed  at  quantifying  Fe:C  uptake 
ratios  as  a  function  of  ambient  [Fer]  that  are  within  the  range  observed  in  the  surface 
ocean  would  help  in  this  regard. 

The  source(s),  sink(s)  and  chemical  characterization  of  Fe-binding  ligands  are  not 
well  known,  resulting  in  the  parameterization  of  ligand  chemistry  with  an  imposed, 
finite  [Lx].  Experimental  work  aimed  at  studying  the  ligand  cycle  would  provide 
information  to  parameterize  ligand  chemistry  in  a  more  mechanistic  manner  than  is 
presently  done. 

There  is  a  pressing  need  for  more  iron  measurements,  especially  in  the  deep  wa- 
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Figure  7-1:  Schematic  diagram  of  ecosystem  model.  Figure  courtesy  of  S.  Dutkiewicz. 

ters.  As  shown  in  Chapter  1,  there  are  only  a  handful  of  measurements  at  depths 
greater  than  2500m,  making  it  very  difficult  to  constrain  the  model  at  depths.  Based 
on  model  results,  I  suggest  measurements  are  needed  in  the  Southern  Hemisphere, 
particularly  in  the  Southern  Ocean  and  the  southern  Indo-Pacific  basin.  My  model 
results  (Chapter  5)  suggest  that  iron  has  a  seasonal  cycle  at  the  surface.  Time-series 
of  surface  measurements  at  both  high  and  low  flux  regions  would  be  useful  in  this 
regard. 

7.1.2  On-going  Work 

The  model  does  a  poor  job  of  predicting  surface  [Fer]  in  high  dust  flux  regions  such 
as  the  North  Atlantic.  I  am  working  to  improve  the  model’s  performance  in  such 
regions  by  including  solubility  limits  of  iron  and  processes  in  surface  waters  that  may 
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serve  as  a  sink  for  iron  in  surface  waters,  such  as  colloid  formation  and  aggregation 
of  iron. 

Replacing  the  simplified  export  parameterization  with  an  explicit  ecosystem  model 
would  distinguish  the  varying  Fe  requirements  of  different  sized  phytoplankton.  In 
addition,  recent  work  has  shown  that  iron  limitation  increases  Si:N  uptake  ratios 
(Hutchins  and  Bruland,  1998;  Takeda,  1998).  Accounting  for  such  a  response  may 
improve  the  model’s  ability  to  predict  modern  surface  [PO4]  and  would  have  impli¬ 
cations  for  an  increased  dust  flux  scenario.  Currently,  Stephanie  Dutkiewicz,  Mick 
Follows  and  I  are  coupling  a  simple  ecosystem  model  with  the  iron  model  presented  in 
this  thesis.  The  ecosystem  model  consists  of  three  nutrient  pools:  nitrogen,  silica,  and 
iron,  as  well  as  two  phytoplankton  size  classes  to  account  for  the  differing  physiologi¬ 
cal  requirements  (Figure  7-1).  The  model  successfully  captures  primary  productivity 
(Figure  7-2).  We  plan  to  continue  developing  these  models  and  perform  sensitivity 
tests  to  explore  implications  of  global  change. 

Given  the  importance  of  the  Southern  Ocean  with  regards  to  climate  change,  more 
work  needs  to  be  done  in  this  area  to  better  understand  its  role  in  absorbing  CO2-  As 
discussed  in  the  preceding  paragraphs,  improving  the  surface  macronutrient  and  iron 
concentrations  during  the  modern  will  help  in  determining  the  Southern  Ocean’s  role 
in  past /future  C02  absorption.  Furthermore,  a  better  understanding  of  dynamics  in 
this  region  is  also  necessary.  Due  to  the  importance  of  diatoms  in  the  Southern  Ocean, 
it  is  also  important  to  couple  iron  and  silica  cycles.  To  better  constrain  the  Southern 
Ocean  results,  I  am  working  with  Taka  Ito,  who  has  developed  a  dynamical  model  of 
the  Southern  Ocean,  to  study  tracer  transport.  We  are  also  coupling  the  ocean  model 
to  a  simple  atmospheric  model  to  make  projections  regarding  pC02  drawdown. 
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Total  Phytoplankton 


Figure  7-2:  Comparison  of  modeled  chlorophyll  (top  panel)  with  chlorophyll  calcu¬ 
lated  from  SeaWiFs.  Figures  courtesy  of  S.  Dutkiewicz. 
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7.2  Final  Comments 


The  purpose  of  this  thesis  was  to  further  the  understanding  of  iron  cycling,  an  element 
known  to  limit  biological  productivity  in  upwelling  regions  of  the  global  ocean.  The 
results  of  this  thesis,  in  agreement  with  other  models,  suggests  that  a  ten  fold  increase 
in  global  dust  flux  could  drawdown  as  much  as  0.5/xM  PO4  in  the  Southern  Ocean, 
approximately  1/3  of  the  excess  [PO4],  clearly  not  the  only  process  affecting  pC02- 
Working  on  modeling  the  iron  cycle  has  made  it  very  clear  to  me  that  the  ocean 
nutrient  cycle  is  extremely  complex  and  there  are  many  open  questions  regarding  our 
understanding  of  the  system.  Trying  to  tinker  with  the  cycling  of  iron  to  mitigate 
increased  pC02  may  lead  to  many  unintended  consequences.  As  Rachel  Carson,  the 
noted  marine  scientist  and  environmentalist  stated,  “the  ‘control  of  nature’  is  a  phrase 
conceived  in  arrogance,  born  of  the  Neanderthal  age  of  biology  and  the  convenience 
of  man.”  Rather,  “the  more  clearly  we  can  focus  our  attention  on  the  wonders  and 
realities  of  the  universe  about  us,  the  less  taste  we  shall  have  for  destruction.” 
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